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ABSTRACT 


Requirements  for  accurate  bend-testing  of  four-point  and  three-point 
beams  of  rectangular  cross-section  are  outlined.  The  so-called  simple 
beam  theory  assumptions  are  examined  to  yield  beam  geometry  ratios  that 
will  result  in  minimum  error  when  utilizing  elasticity  theory.  Factors 
that  give  rise  to  additional  errors  when  determining  bend  strength  are  ex¬ 
amined,  such  as:  wedging  stress,  contact  stress,  load  mislocation,  beam 
twisting,  friction  at  beam  contact  points,  contact  point  tangency  shift, 
and  neglect  of  corner  radii  or  chamfer  in  the  stress  determination. 

Also  included  are  the  appropriate  Weibull  strength  relationships  and  an  ! 
estimate  of  errors  in  the  determination  of  the  Weibull  parameters  based  on 
sample  size.  Error  tables,  based  on  many  of  the  previously  mentioned  fac¬ 
tors,  are  presented. 

Such  data  are  utilized  to  rationally  arrive  at  a  1/3-four-point  and  a 
three-point  beam  of  a  specific  geometry  ratio  and  dimension  that  will  result 
in  minimization  of  errors  during  flexure  testing.  Although  such  beams  are 
recommended  for  a  military  standard,  they  are  termed  reference  standard 
beams.  Beams  of  other  loading  configurations  and  dimensions  are  considered 
acceptable  if  statistical  correlation  to  one  of  the  reference  standards  is 
determined,  and  if  error  spans  are  reported  when  presenting  flexure  test 
data, 
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INTRODUCTION 


There  has  been  an  increase  in  interest  and  activity  in  recent  years  in  both 
the  research  and  development  of  ceramic  materials  and  their  practical  application 
to  engineering  structures.  Much  of  the  impetus  is  due  to  the  Defense  Advanced 
Research  Projects  Agency  (DARPA)  program1  where  the  objective  is  to  develop  ceramic 
turbine  engines  for  vehicular  use  and  electric  power  generation.  It  was  recognized* 
during  the  course  of  this  program  that  technical  progress  would  be  improved  by  having 
available  unified  standard  testing  methods  so  that  the  interchange  of  test  data  be¬ 
tween  involved  organizations  would  have  the  same  reference  base.  This  requirement 
is  vital  because  ceramic  materials  need  particular  kinds  of  mechanical  tests  suit¬ 
able  to  their  brittle  nature  so  that  valid  test  data  will  result. 

Action  was  initiated  to  determine  what  kind  of  standards,  if  indeed  any,  were 
available  for  mechanical  testing  of  brittle  materials. +  It  was  determined  from  this 
study  that  "of  some  80-odd  existing  ceramic  standards  reviewed  (c.f.  NBS  SP  329,  and 
Index  of  U.S.  Voluntary  Engineering  Standards  and  DoD  Index  of  Specifications  and 
Standards),  only  one  was  found  suitable  for  testing  structural  ceramics.  An  in-house 
review  and  survey  of  cognizant  government  and  nongovernment  activities  confirmed  both 
the  lack  of  test  method  standards,  and  the  real  need  for  such  standards." 

As  a  result  of  the  above  exercise,  a  tentative  unapproved  set  of  standards^  was 
prepared  by  the  Army  Materials  and  Mechanics  Research  Center  (AMMRC)  and  distributed 
to  interested  and  involved  organizations.  This  set  of  unofficial' standards ,  which 
included  such  test  methods  as  flexure,  tension,  creep,  stress  rupture,  fatigue,  and 
spin  testing,  was  discussed  at  several  meetings  of  government  and  industry  repre¬ 
sentatives.  A  number  of  worthwhile  suggestions  evolved.  However,  it  was  apparent 
that  these  tentative  standards  could  be  improved,  and  such  comments  were  invited. 

These  unofficial  standards  have  remained  unchanged  and  unapproved.  Recently, 
however,  interest  was  revived  at  AMMRC  in  finalizing  standard  tests  for  brittle 
materials.  It  was  viewed  that  the  original  tentative  standards,  dated  2  April  1973, 
represented  the  ideal  goal  but  were  far  too  inclusive  to  realistically  establish 
testing  requirements  which  would  provide  valid  results  at  this  time.  It  had  also 
been  decided  that  each  subject  should  be  separately  treated  as  sufficient  test 
methodology  knowledge  evolved.  Since  flexure  testing  to  determine  bend  strength 
met  this  criterion,  it  was  decided  to  standardize  such  a  test.  Other  specific  recom¬ 
mendations  regarding  flexure  testing  resulting  from  the  AMMRC  committee  were  proposed 
and  are  as  follows: 

The  standard  at  the  present  time  would  consider  only  beam  testing  methodology. 

The  test  would  have  to  conform  within  the  requirements  of  simple  beam  theory, 

thus  the  test  would  be  restricted  to  a  temperature  environment  of  70°F  or  less. 


*The  Ford  Motor  Company,  in  a  letter  dated  31  April  1972,  cited  the  problem  of  nonexisting  standards  for  testing  structural  ceramic 
materials  and  appealed  to  AMMRC  to  initiate  action  to  develop  appropriate  standards. 
fTrip  to  Ford  Motor  Company,  Dearborn,  Michigan,  on  9  May  1973,  by  H.  F.  Campbell,  S.  Acquaviva,  and  L.  G.  MacDonald, 
dated  15  May  1973. 

♦Army  Materials  and  Mechanics  Research  Center,  "Military  Standards  -  Test  Methods  for  Structural  Standards,”  2  April  1973. 

1.  VAN  REUTH,  E.  C.  The  Advanced  Research  Projects  Agency’s  Gas  Turbine  Program.  Proc.  of  the  Second  Army  Materials 
Technology  Conference  -  Ceramics  for  High  Performance  Applications,  J.  i.  Burke,  A.  E.  Gorum,  and  R.  N.  Katz,  ed.. 

Brook  Hill  Publishing  Company,  Chestnut  Hill,  Massachusetts,  1974. 
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Although  two  beam  testing  systems  are  to  be  used,  i.e.,  size,  dimensions, 
loadings,  etc.,  and  called  reference  standards,  other  systems  will  be  accept¬ 
able  as  long  as  strength  relationships  are  determined  and  an  estimate  of  the 
errors  are  reported  along  with  the  test  data  results. 

The  author,  who  had  presented  a  critique*  of  the  tentative  standard,  was  charged 
with  the  task  of  providing  backup  material  in  the  form  of  a  formal  report.  This  work 
is  a  result  of  that  request. 

OBJECTIVE 

The  objective  of  this  report  is  to  recommend  beam  test  systems  such  that  accu¬ 
rate  fracture  strength  utilizing  simple  beam  theory  will  result.  Such  recommendations 
will  be  used  as  the  basis  of  a  military  standard  for  flexure  testing  of  brittle 
materials  within  the  elastic  regime. 


SIMPLE  BEAM  THEORY  ASSUMPTIONS 

The  rectangular  beam  configuration  is  attractive  as  a  strength-test  vehicle 
because  of  its  simple  shape  and  apparent  ease  of  load  application,  as  well  as  analysis 
and  reduction  of  data.  Rods  of  circular  cross  section  are  also  used  in  beam  tests, 
but  usually  for  specialized  testing.  Because  a  beam  of  circular  cross  section  is  not 
as  frequently  used  as  the  rectangular  beam,  the  circular  cross  section  is  not  consi¬ 
dered  a  reference  standard  candidate;  therefore  only  the  rectangular  cross  section  is 
examined  in  the  discussions  to  follow. 

A  critical  review  of  simple  beam  theory  assumptions  will  yield  ranges  of  geometry 
ratios  by  which  the  theory  can  be  validly  applied.  These  assumptions  are  listed  below, 
as  well  as  their  associated  inferences  in  terms  of  an  error  analysis: 

1.  Transverse  planes  perpendicular  to  the  longitudinal  axis  of  the  beam 
remain  plane  after  the  beam  is  bent. 

2.  The  modulus  of  elasticity  in  tension  is  equal  to  the  modulus  of  elasticity 
in  compression.  Also,  the  beam  material  is  isotropic  and  homogeneous. 

3.  The  maximum  deflection  must  be  small  compared  to  the  beam  depth. 

4.  The  beam  must  deflect  normally  under  elastic  bending  stresses  but  not 
through  any  local  collapse  or  twisting. 

Each  of  the  above  assumptions  is  examined  in  detail,  where  possible,  so  that  the 
required  rectangular  beam  geometry  ratios  can  be  determined  as  a  function  of  the 
associated  errors. 

Assumptions  1  and  2  together  imply  that  stress  and  strain  are  proportional  to 
the  distance  from  the  neutral  axis,  and  the  stress  does  not  exceed  the  proportional 
limit  of  the  material.  These  assumptions  disregard  the  effect  of  any  shearing  resist¬ 
ance  and  make  impossible  the  use  of  the  flexure  formula  for  curved  beams  of  large 
curvature. 


•Memorandum  •  Review  of  "Military  Standard  Test  Method  for  Structural  Ceramics,”  by  F.  I.  Baratta,  AMMRC,  December  1973, 
revised  February  1974. 
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Assumption  1  and  the  above  implication  infer  that  there  is  a  linear  relationship 
of  stress  between  the  fibers  of  the  beam.  Timoshenko  and  Goodier2  give  the  more  exact 
distribution  of  stress  as  the  function  of  the  distance  from  the  neutral  axis  for  a 
simply  supported  beam  loaded  by  distributed  pressure.  Although  the  candidate  beams 
to  be  considered  here  will  be  four-point  and  three-point  loaded  beams,  these  different 
loading  conditions  should  not  yield  markedly  different  results  from  those  taken  from 
the  above  reference.*  From  Reference  2  an  estimate  of  the  error  incurred  when  assuming 
a  linear  stress  distribution  can  be  obtained  from  the  following  equation: 


The  above  equation  is  approximately  applicable  for  both  four-point  and  three-point 
loading.  The  beam  geometries  and  loading  configurations  considered  are  shown  in 
Figures  la  and  2a.  The  error  resulting  from  assuming  a  linear  stress  distribution 
when  in  fact  nonlinearity  exists  is  given  in  Table  1.  (Tables  begin  on  page  25.) 

Regarding  the  assumption  that  the  modulus  of  elasticity  in  tension  is  equal  to 
that  in  compression,  Ej  =  Ec,  Chamlis3  has  derived  in  closed  form  the  solution  for 
the  tensile  bending  stress  when  E-p  f  Ec.  After  some  manipulation  of  appropriate  formu¬ 
las  in  Reference  3,  the  tensile  stress  due  to  bending  is  given  by: 

ox  =  (ob/2)[l  +  (Et/Ec)2],  (2) 

for  both  the  four  point  and  the  three-point  loaded  beams.  The  resulting  percent  error 
is  given  in  Table  2. 

Although  the  errors  associated  with  neglecting  to  account  for  anisotropy  and 
nonhomogeneity  of  the  test  material  are  not  considered  here,  they  arc  briefly  mentioned 
in  the  following  paragraphs  so  that  the  reader  will  be  aware  of  such  possibilities. 

If  the  beam  is  anisotropic,  the  bending  stress  formula  is  exactly  the  same  as 
the  elementary  theory  except  that  the  application  of  a  bending  moment  produces  a  twist¬ 
ing  moment,  and  vice  versa.  According  to  Lekhnitskii ,  '*  determining  the  accompanying 
shear  stress  produced  by  bending  a  rod  of  rectangular  cross  section,  having  only  one 
plane  of  elastic  symmetry  normal  to  the  axis,  is  very  complicated.  (Composite  and 
crystal  structures  are  excluded  here  as  test  materials.)  If  the  degree  of  anisotropy 
for  ceramic  material  is  slight,  it  may  be  permissible  to  assume  that  the  error  when 
ignoring  this  effect  on  the  fracture  stress  will  also  be  small. 

Nonhomogeneity  of  the  test  material  infers  variation  of  the  elastic  modulus. 

It  has  been  observed!-  that  in  plates  of  hot -pressed  silicon  nitride,  the  modulus  of 
elasticity  at  the  surface  is  several  percent  different  than  that  of  the  center.  This 


1  *  (4/l5)(-C7 3-)"’ 


•Recent  finite  element  computer  results  by  H.  I  aou  and  R.  J.  fl.  bollard  of  the  University  of  Washington  indeed  show  this  to  be  true. 
See  flexure  Test  Method  T10,  Proceedings  of  the  Sixth  Army  Materials  Technology  Conference.  “Ceramics  for  High  Performance 
Applications  -  III  -  Reliability,"  held  at  Orcas  Island.  Washington.  July  1979. 
fPrivate  discussion  with  P.  M.  Lcnoe,  AMMRC. 

2.  TIMOSHENKO,  S..  and  GOODII  R,  J.  N.  Theory  of  Elasticity.  2nd  I'd..  McGraw-Hill  Book  Co.,  Inc..  New  York.  1951. 

3.  CHAMLIS,  C.  C.  Analysis  of  ThrcePoint-Bend  Test  for  Materials  with  Unequal  Tension  and  Compressive  Properties.  NASA  TN 
D7572.  March  1974. 

4.  I.IKHNITSKII.  S.  G.  Theory  of  Elasticity  of  an  Anisotropic  Elastic  Body.  Holdcn-Day  Scries  in  Mathematical  Physics, 

J.  L.  Brandstatter,  cd„  1963,  p.  204. 
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is  also  an  area  in  which  further  analysis  will  be  required  to  assess  the  error  appli¬ 
cable  to  four-point  and  three-point  loaded  beams  when  the  modulus  of  elasticity  varies 
through  the  beam  thickness. 

If  a  rectangular  beam  .has  initial  curvature  pc,  the  error  can  be  determined  from 
an  analysis  provided  by  Timoshenko.5  The  general  bending  stress  ox  in  a  curved  beam 
due  to  a  pure  moment  is  given  by  the  following: 

ax  =  ac(Mb/bdpc)  (3) 

where 

(d/2pc)  -  (ec/pc) 

ac  =  -  (3a) 

(ec/pc)[l  -  (d/ 2pc)  ] 

ec/pc  =  [(d/pc)7l2][l  +  (d/pc)7l5].  (3b) 

Since  the  bending  stress,  according  to  simple  beam  theory,  is  ob  =  6Mb/bd2,  and 
putting  ob  in  the  same  terms  as  (3)  above,  we  have: 

°b  =  %(Mb/bdPc7  (7 

where 

%  =  6(pc/d).  (4a) 

The  percent  error  J  for  a  bent  beam  of  rectangular  cross  section  and  of  beam- 
to-depth  initial  curvature  Pc/d  resulting  in  a  neutral  axis  shift  of  ec/pc  is: 

e  =  100 [  (ab  -  ac)/ac] .  (5) 

The  resulting  error  for  a  beam  of  rectangular  cross  section  bent  by  a  pure  moment  as 
obtained  from  (5)  is  given  in  Table  3  as  a  function  of  initial  curvature.  It  is  as¬ 
sumed  that  an  analogous  analysis  applied  to  a  three-point  loaded  beam  would  produce 
similar  results. 

The  validity  of  the  inference  that  the  strain  is  proportional  to  the  distance 
from  the  neutral  axis  is  dependent  upon  the  ratio  of  the  beam  width  to  its  depth. 
Anticlastic  curvature  of  rectangular  beams  or  plates  with  intermediate  ratios  of  b/d 
can  lead  to  erroneous  results  using  simple  beam  theory;  see  Timoshenko.6  Of  course, 
if  the  beam  can  be  considered  infinite  in  width,  like  a  plate,  the  correction  of  the 
bending  stress  is  simply  1/(1  -  v2).  The  question  arises  as  to  what  ratios  of  b/d 
are  appropriate  for  the  application  of  simple  beam  theory.  Ashwell7  examined  in  detail 
the  anticlastic  curvature  of  rectangular  beams  and  plates  and  provided  the  answer  to 
this  question.  The  pertinent  formulas  taken  from  Reference  7  are  given  in  Appendix  A. 
These  equations  were  applied  to  ceramic  materials  with  Poisson's  ratio  v  equal  to  0.25 
and  the  ratio  of  Young's  modulus  to  fracture  stress  E/ob  of  1  *  103  to  determine  the 
percent  error*  using  simple  beam  theory  as  a  function  of  b/d  which  is  shown  in  Table  4. 


•Ashwell  considered  a  beam  bent  by  a  constant  moment  analagous  to  the  four-point  beam  loading  case,  which  should  represent  a 
conservative  bound  on  b/d  for  the  three-point  beam,  as  well. 

5.  TIMOSHENKO,  S.  Strength  of  Materials.  Parts  I  and  II,  D.  Van  Nostrand  Co.,  Inc.,  New  Jersey,  1976. 

6.  TIMOSHENKO,  S.  Letter  to  the  Editor.  Mechanical  Engineering,  v.  45,  no.  4,  April  1923,  p.  259-260. 

7.  ASHWELL,  D.  G.  The  Anticlastic  Curvature  of  Rectangular  Beams  and  Plates.  J.  Roy.  Aero.  Soc. ,  v.  54,  1950,  p.  708-715. 
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If  the  maximum  deflection  is  not  small  compared  to  the  beam  depth,  linear  beam 
theory  cannot  be  employed  without  an  error.  West8  examined  large  deflections  of  three- 
point  loaded  beams,  and  from  such  results  a  definitive  ratio  of  beam  length  to  depth 
can  be  determined  for  valid  application  of  simple  beam  formulas.  Since  for  most  brittle 
materials  values  of  E/ob  range  from  approximately  0.5  x  103  to  1  x  103,the  former  val¬ 
ue  was  used  to  compute  the  percent  error  because  it  would  yield  the  largest  error.  Al¬ 
though  the  analysis  was  applied  to  a  three-point  loaded  beam,  the  method  was  extended 
to  determine  errors  for  four-point  loaded  beams  as  well.  The  results  of  the  calcula¬ 
tions  using  the  mentioned  analysis®  are  presented  in  Table  5,  which  gives  errors  for  a 
1/3-four-point,  a  1/4-four-point,  and  a  three-point  loaded  beam  as  a  function  of  L/d. 

It  is  implicit  in  the  assumptions  given  in  Reference  8  that  the  loads  and  moments 
are  applied  to  the  beam  in  an  ideal  manner  with  no  friction  occurring  between  the  load 
application  points  and  the  beam.  Ritter  and  Wilson9  have  determined  a  beam  length-to- 
depth  limit  based  on  the  minimization  of  friction  effects  when  large  deflections  occur. 
The  friction  effect  considered  is  that  which  gives  rise  to  a  moment  caused  by  the  slope 
at  the  load  application  point.  Not  considered  in  the  analysis9  are  the  effects  of  fric¬ 
tion  due  to  a  moment  acting  out  of  the  neutral  plane  of  the  beam,  lateral  contraction 
or  extension,  and  changes  in  moment  arms  due  to  contact  point  tangency  shift.  These 
factors  will  be  discussed  later. 

Returning  to  the  results  of  Reference  9,  an  inequality  for  the  four-point- loaded 
beam  which  provides  a  limit  is  given  in  the  following: 

(L/d  -  a/d)/(E/ob)  <  0.3  (6) 

to  insure  negligible  nonlinear  deflections  and  friction  effects.  The  value  of  0.3  was 
obtained  from  limiting  the  slope  to  less  than  15°  between  the  beam  in  the  loaded  and 
unloaded  positions  at  the  outermost  support  point.  If  the  minimum  value  of  E/ob  is 
chosen  to  be  0.5  x  10 3 ,  then  we  determine  that  for  a  four-point  loaded  beam  (6)  becomes: 

L/d  -  a/d  £  150.  (7) 

It  is  noted  from  Table  5  that  neglecting  beam  deflections  resulted  in  greater 
error  in  calculation  of  bending  stress  for  the  four-point  loaded  beam  than  for  the 
three-point  loaded  beam.  For  conservatism,  therefore,  it  will  be  assumed  that  (7)  is 
applicable  to  the  three-point  loaded  beam  as  well,  with  a/d  =  0.  Thus  (7)  becomes 

L/d  £  150.  (8) 

It  appears  that  these  limits  are  compatible  with  those  values  given  in  Table  5  such 
that  reasonable  L/d  ratios  can  be  chosen  that  will  result  in  small  errors  when  mini¬ 
mizing  deflection. 

One  of  the  last  requirements,  no  buckling  of  the  beam,  is  easily  fulfilled  for 
ceramic  materials  with  beam  dimensions  of  practical  test  configurations.  The  reader 
can  readily  verify  this  statement  by  referring  to  Timoshenko  and  Gere.10 

Accuracy,  which  is  inferred  in  the  above  restrictions,  is  also  dependent  upon  the 
manner  of  load  application,  beam  geometry,  loading  fixtures,  and  surface  preparation. 
Although  specimen  size  will  not  affect  accuracy  except  for  extremely  small  geometries, 

8.  WEST,  D.  C.  Flexure  Testing  of  Plastics.  Exp.  Mech.,  v.  21,  no.  2,  July  1964. 

9.  RITTER,  J.  E„  and  WILSON,  W.  R.  D.  Friction  Effects  in  Four-Point  Bending.  ASLE  Transactions,  v.  18,  no.2,  p.  130-134, 
presented  at  the  29th  annual  meeting,  April  28-May  2,  1974. 

10.  TIMOSHENKO,  S.,  and  GERE,  J.  M.  Theory  of  Elastic  Stability.  McGraw-Hill  Book  Co.,  Inc.,  New  York,  1961. 
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it  will  alter  the  magnitude  of  the  stress  level  at  failure,  and  this  must  also  be  con¬ 
sidered.  These  subjects  are  discussed  in  the  following  paragraphs,  and  guidelines  for 
specimen  geometry  and  minimization  of  errors  are  provided. 

First  to  be  considered,  however,  are  the  merits  of  a  four-point  beam  loading 
system  as  compared  to  the  three-point  beam  loading  system. 


FOUR-POINT  AND  THREE-POINT  LOADING 

The  bending  moment,  from  which  the  desired  fracture  stress  is  computed  in  an 
idealized  four-point  beam  loading  system,  as  shown  in  Figure  la,  is  constant,  and 
there  are  no  horizontal  or  vertical  shear  stresses  within  the  inner  span.  However, 
the  bending  moment  in  an  idealized  three-point  beam  loading  system,  shown  in  Figure  2a, 
is  linearly  dependent  upon  the  distance  from  the  nearest  support  to  the  fracture  ori¬ 
gin,  and  thus  requires  an  additional  distance  measurement  to  determine  the  fracture 
stress.  Also,  the  shear  stresses  for  the  three-point  beam  loading  system  are  devel¬ 
oped  over  the  full  span,  thus  deviating  from  the  ideally  sought  uniaxial  stress  state 
present  in  the  four-point  beam  loading  system. 

Wedging  stresses5  occur  under  all  points  of  load  application  during  flexure 
testing  of  beams.  The  effect  of  the  wedging  stress  occurring  at  the  inner  load  points 
of  a  four-point  beam  test  is  to  cause  a  deviation  from  the  idealized  calculated  constant 
stress  at  the  two  local  regions.  However,  if  the  ratio  of  half  the  distance  between  the 
outer  span  L  and  inner  span  i,  called  a,  to  beam  depth  d  is  great  enough,*  the  stress 
reduction  will  not  only  be  small  but  will  decay  rapidly,  and  the  stress  predicted  by 
simple  beam  theory  will  be  developed.  Yet,  the  maximum  stress  computed  by  simple  beam 
formula  for  the  three-point  beam  system  is  never  attained.  The  actual  maximum  stress 
occurs  at  a  short  distance  either  side  of  the  center  of  the  load  application  point, 
which  can  cause  fracture  at  these  sites,  rather  than  at  the  center,  according  to 
Rudnick  et  al.11  This  observation  has  also  been  confirmed  by  Oh  and  Finnic,1 ?  where 
only  for  a  material  with  no  scatter  in  strength  will  the  fracture  location  of  a  three- 
point  loaded  beam  be  theoretically!  located  at  the  central  load  point. 

Brittle  materials  are  affected  by  size.  Compensation  can  be  realized  through  the 
use  of  statistical  analysis  offered  by  Weibull.13  Although  the  four-point  beam  system 
assures  a  simple  stress  state  which  is  easier  to  analyze1  than  the  more  complex  biaxial 
stress  state  associated  with  the  three-point  beam  specimen,  this  will  be  less  of  a  con¬ 
sideration  if  the  beam  is  designed  properly.  Nevertheless,  the  three-point  loaded  beam 
system  is  preferred  when  investigating  material  or  process  development,  because  of 
smaller  specimen  size,  or  when  attempting  to  pinpoint  fracture  origin  location. =F  On 
the  other  hand,  the  four-point  loaded  beam  is  preferred  when  determination  of  strength 
for  design  purposes  is  desired,  because  the  center  span  is  uniaxially  stressed,  i.e., 
no  shear  stresses  exist.  It  is  concluded  that  each  of  these  systems  is  suited  for  a 
particular  application  and  each  has  different  advantages  and  disadvantages.  Thus  both 
types  of  loading  systems,  the  four-point  beam  system  shown  in  Figure  1  and  the  three- 
point  beam  system  shown  in  Figure  2,  will  be  considered  as  reference  standards. 

•This  requirement  will  be  discussed  subsequently. 

tin  Reference  12,  the  authors  considered  only  a  statistical  analysis  and  ignored  wedging  stress  considerations. 
fPrivate  communication  with  R.  W.  Rice  of  N.  R.  L. 

11.  RUDNICK,  H„  MARSCHALL,  C.  W„  DUCKWORTH,  W.  H.,  and  FNR1CK ,  B.  R.  The  Evaluation  and  Interpretation  of  Mechanical 
Properties  of  Brittle  Materials.  Al'ME  TR  67-316,  April  1968. 

12.  OH,  11.  L..  and  FINNIK,  1.  On  the  Location  of  Fracture  In  Brittle  Solids  ■  I.  Due  to  Static  Loading  lnt.  J.  of  Fracture  Mechanics, 
v.  6,  no.  3,  September  1970,  p.  287-300. 

13.  WEIBULL,  W.  Statistical  Theory  of  Strength  of  Materials.  Royal  Swedish  Institute  for  Engineering  Res.,  Proc.  no.  151,  1939,  p.  1-45. 
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Each  of  these  beam  systems  will  be  subjected  to  external  influences  which  will 
affect  the  accuracy  of  the  test  results.  These  external  influences,  directly  or  indi¬ 
rectly  caused  by  the  application  of  loads  through  the  test  fixtures,  will  lead  to 
either  configuration  constraints  or  errors. 

EXTERNAL  INFLUENCES 

The  major  influence  on  the  accurate  determination  of  flexure  strength  of  a  beam 
in  bending  arises  from  the  application  of  load  through  the  fixtures  to  the  specimen. 

The  idealizations  indicated  in  Figures  la  and  2a  are  rarely  met,  and  usually  tests 
are  conducted  using  a  convenient  rigid  loading  head  and  support  member  as  depicted  in 
Figures  lb  and  2b.  The  constraints  on  either  the  loading  fixture  or  the  specimen 
and/or  errors  resulting  from  such  fixture  designs  are  many.  Such  constraints  or  errors, 
which  are  discussed  in  turn,  are  caused  by: 

1.  load  mislocation 

2.  beam  twisting 

3.  friction 

4.  local  stresses 

5.  contact  point  tangency  shift 

6.  surface  preparation 

1.  Load  Mislocation 

a.  Four-point  loaded  beams 

When  calculating  bending  stress  by  simple  beam  theory  formula  for  four-point 
loaded  beams,  it  is  usual  to  assume  that  the  moment  within  the  inner  span  l  is  constant. 
However,  if  a  loading  head  that  can  only  translate  is  used,  as  idealized  in  Figure  lb, 
it  is  impossible  to  attain  this  idealized  moment  condition  when  Xj  /  X3  -  x2;1114  this 
is  shown  in  Figure  lc.  The  ratio  of  ax/ob,  from  Appendix  B,  is: 

ax/ab  =  - — -  xx/a  .  (9) 

.  (Pz  +  Ps)/2 

The  loads  and  distances  are  also  shown  in  Figure  lc ,  and  a  is  the  value  of  ai  with 
perfect  load  location.  The  error  is  magnified  by  the  ratio  of  xi/a.  (Of  course,  if 
Pi  =  P2  =  P3,  which  implies  exact  location  of  the  points  of  load  application,  there  is 
no  error.)  In  order  to  estimate  the  magnitude  of  such  an  error  it  was  assumed  in 
Appendix  B  that  the  upper  two  load  points  in  Figure  lc  were  at  a  fixed  distance  x2-xi  =  f. 
and  were  constrained  to  translate  vertically  during  loading,  and  that  the  loading  head 
would  be  located  such  that  Xi  /  X3  -  x2.  This  method  of  loading,  being  the  most  conveni¬ 
ent,  is  usually  adopted  by  many  investigators,  and  therefore  the  resulting  error  deter¬ 
mination  is  not  unrealistic. 

The  analysis  was  accomplished  by  simply  enforcing  the  condition  that  the  displace¬ 
ment  at  X!  must  be  equal  to  the  displacement  at  x2  in  the  deflection  equation.  This 
results  in  the  following  relationships  between  ox  and  ab  in  terms  of  the  load  eccen¬ 
tricity  ratio  e/L: 

14.  HOAGLAND,  R.  G.,  MARSCHALL,  C.  W.,  and  DUCKWORTH,  W.  H.  Reduction  of  Krrors  in  Ceramic  Bend  Tests.  J.  Amer. 

Cer.  Soc.,  v.  59.  no.  5-6,  May-June  1976,  p.  189-192. 
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°x/°b= 


f 


[(e/L  +  a/L) / (a/L) ]  [  1- (e/L  +  a/L)  -  £/L]  { (Jl/L)  [2- (e/L  ♦  a/L)]-2[l-(e/L  ♦  a/L)]2} 


3 (e/L  +  a/L) [1-  Jl/L  -(e/L+  a/L)]-(l-  Jl/L)2 
where  the  parameters  a,  i  and  L  are  shown  in  Figure  1. 


(10) 


Most  workers  in  the  testing  field  utilize  either  a  1/3-point  (a/L  =  1/3  and 
1/ L  =  1/3)  or  a  1/4-point  (a/L  =  1/4  and  Jl/L  =  1/2)  loading.  Thus  by  substitution  of 
these  parameters  into  Equation  10,  we  obtain: 


(°x/°b) 


*/'L=l 


[3(e/L)+l] (1/3  -  e/L)[l/3(5/3  -  e/L)-2(2/3  -  e/L)2] 
[3(e/L)+l] (1/3  -  e/L)  -  4/9 


(ID 


'(°X/0b^/L=l 


[4 (e/L) +1] (1/4  -  e/L) [l/2(7/4  -  e/L)-2(3/4  -  e/L)2] 
[3(e/L)  +  3/4] (1/4  -  e/L)  -  1/4 


(12) 


The  reader  is  cautioned  that  for  given  values  of  Jl/L  there  exists  a  limit  on  e/L 
in  (10),  (11),  and  (12);  that  is,  ai  can  be  such  that  either  P2  or  P3  =  0  because  the 
test  system  changes  from  four-point  to  an  eccentric  three-point  loading.  (See  App.  B.) 

The  error,  defined  as  [(ajj-ox) /ax] 100,  was  determined  from  (11)  and  (12)  for  the 
1/3-point  and  1/4-point  loaded  beams  and  is  shown  in  Tables  6  and  7  as  a  function  of 
e/L.  Even  though  these  tables  show  ±e/L  values,  only  negative  values  of  e/L  were  con¬ 
sidered  in  (11)  and  (12)  because  when  e/L  <0,  ox  >  %.  Tables  6  and  7  show  that  for 
corresponding  e/L,  when  a.x/1  t  a2/L,  the  1/3-point  loading  system  results  in  lesser 
error  than  the  1/4-point  loading  system.  Also,  in  accordance  with  the  above  discus¬ 
sion,  e/L  in  Tables  6  and  7  is  limited  to  a  range  of  ±0.0443  and  ±0.0465.  The  errors 
indicated  in  these  tables  can  be  minimized  by  designing  the  loading  fixture  so  that 
the  inner  and  outer  spans  are  independently  fixed.  Also,  the  inner  span  should  be 
designed  with  accurate  location  adjustment  and  allowed  to  pivot  as  recommended  by 
Hoagland  et  al.1H 

b.  Three-point  loaded  beams 

The  same  type  of  analysis  was  applied  to  the  three-point  loaded  beam.  The 
ideal  system  is  shown  in  Figure  2a  and  the  system  analyzed  is  shown  in  Figure  2b.  The 
resulting  percent  error  due  to  an  eccentric  load  application  is  given  by 


e  =  (1/4)  ( - - - )  (  — —  V  100. 

\  1-ai/L  )  \ai/L  /  J 

The  percent  error  as  a  function  of  e/L  is  given  in  Table  8. 

Notice  that  the  percent  errors  given  in  Table  8  are  always  positive,  and  when 
the  load  application  point  is  misplaced,  such  errors  are  much  less  than  those  of 


equivalent  e/L  values  shown  in  Tables  6  and  7  for  the  four-point  loaded  beams. 

Notice  also  that  when  ±e/L  =  0.500,  the  error  is  infinite,  i.e.,  the  three-point 
loading  model  is  no  longer  valid. 

2.  Beam  Twisting 

A  net  torque  can  result  from  line  loads  being  nonuniform  or  nonparallel  between 
pairs  of  load  contact  points  or  if  the  cross  section  of  the  specimen  is  skewed  over 
its  length.11'14 

Such  a  condition  is  shown  schematically  in  Figure  3  for  a  four-point  bending 
specimen.  The  error  due  to  twisting  has  been  estimated14  for  both  the  plane  stress 
and  plane  strain  conditions  by  examining  the  maximum  principal  stress  due  to  bending 
and  torsion  and  comparing  it  to  the  bending  stress.  Only  the  details  of  the  plane 
stress  analysis  were  given  in  Reference  14,  and  bottoming  of  the  specimen  on  the 
fixture  was  not  considered.  For  the  sake  of  completeness,  the  plane  strain  condi¬ 
tion  and  bottoming  are  considered  in  the  analysis  given  in  Appendix  C. 

The  maximum  principal  stress  for  either  a  skewed  four-point  or  three-point 
beam  in  bending,  considering  a  plane  strain  condition,  is  given  by: 

°n  =  (V2){1+v+(1/3k2)[(n/ll,/b)2  +  9k22(l-v)2]2)  (14a) 

max 

where  is  the  apparent  bend  strength  and  a •  is  either  equal  to  "a”  for  four-point 
bending  or  equal  to  L/2  for  three-point  bending.  Also: 

n  =  [SkiCE/ahJ/CUvmCd/LrHs  +  (d/A')*F](b/l')  (14b) 

where  for  Case  I:  n  =  1,  failure  occurs  prior  to  bottoming  of  the  specimen  in  the 
loading  fixture,  and  for  Case  II:  n  <  1,  failure  occurs  after  bottoming. 

The  factors  k!  and  k2,  obtained  from  Reference  2  and  given  in  Table  9,  are 
numerical  values  associated  with  the  torsional  stress  component  which  are  dependent 
on  the  ratio  of  b/d.  The  measured  angle  of  twist  (or  skew  angle)  along  the  total 
length  Ly  of  the  specimen  is  <|>s  (see  Figure  3)  .and  along  the  fixture  from  one  support 
point  to  the  adjacent  load  point,  it  is  <j>p . 

The  maximum  principal  stress*  as  given  by  (14a)  can  be  utilized  to  determine 
the  percent  error  for  various  ratios  of  n,  2’/b,  and  b/d.  This  was  accomplished  and 
is  shown  in  Table  10.  Notice  that  the  range  of  n  varies  from  0.20  to  1.00.  It  is 
expected  that  if  bottoming  does  not  occur  prior  to  fracture  because  of  an  excessive 
twist  angle,  the  maximum  ratio  of  n  that  can  be  attained  is  1.0  and  thus  the  tables 
do  not  accommodate  n  >  1.0.  The  use  of  such  tables  is  demonstrated  in  "Discussion 
of  Errors." 

3.  Friction 

It  has  already  been  shown  in  Table  5,  for  the  two  beam  systems  considered,  that 
the  error  due  to  deflection  will  be  equal  to  or  less  than  if  L/d  <.  25.  It  appears 
that  these  limits  are  well  within  the  geometry  ratios  required  for  the  reference  stan¬ 
dards  to  be  proposed  here.  Therefore,  friction  at  the  load  and  support  points  will 

•The  maximum  principal  str era,  assuming  a  plane  stress  condition,  can  be  realized  by  simply  allowing  Poisson's  ratio  v  to  be  zero  in 
Equation  14a.  However,  this  will  result  in  a  lesser  error  than  when  assuming  the  plane-strain  condition. 
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p  p 


p  p 
(b)  End  View 


Figure  3.  Twisting  of  a  four-point  beam  specimen. 


have  a  negligible  effect  with  regard  to  the  use  of  simple  beam  theory.  This  also 
implies  that  there  will  be  no  effect  from  friction  on  the  contact  tangency  shift. 
(These  factors  will  be  discussed  subsequently.)  However,  friction  will  cause  a 
moment  acting  out  of  the  plane  of  the  beam  that  can  not  be  ignored.  This  factor  is 
considered  in  the  following. 


When  determining  bend  strength  by  simple  beam  theory,  it  is  usual  to  assume 
that  the  supports  and  load  points  are  frictionless,  whereas  in  fact  they  are  not. 

The  presence  of  friction  in  flexure  tests  with  fixed  load  and  support  points  gives 
rise  to  couples  at  such  locations  as  well  as  axial  forces  at  the  neutral  axis  of  the 
beam.  The  net  axial  force  is  relatively  small  and  therefore  is  ignored  here.  How¬ 
ever,  if  the  moment  is  not  corrected  to  account  for  the  couple  in  the  determination 
of  flexure  stress,  an  error  will  result.  Error  equations  adapted  from  the  results* 
available  in  the  literature14-16  are  given  below  for  the  four-point  and  three-point 
loading  systems: 


(15) 


(lb) 


Such  errors  as  defined  by  the  above  equations  can  be  significant,  according  to 
References  14,  16,  and  17.  Newnham16  and  Weil17  reported  that  the  experimental  dif¬ 
ference  in  failure  stress  using  rigid  knife  edges  as  compared  to  roller-type  contact 
points  was  as  high  as  12%  for  silicon  nitride  and  13%  for  graphite. 


•Beam  width  constraint  occurs  also  because  of  friction  transverse  to  the  beam's  long  axis.  However,  this  effect  (see  Newnham'^)  is 
small  and  thus  not  considered  here. 

15.  DUCKWORTH,  W.  H..  et  al.  Mechanical  Property  Tests  on  Ceramic  Bodies.  WADC  TR  52-67,  March  1952,  p.  67-70. 

16.  NEWNHAM,  R.  C.  Strength  Tests  for  Brittle  Materials.  Proc.  of  the  British  Cer.  Soc.,  no.  25,  May  1975,  p.  281-293. 

17.  WEIL,  N.  A.  Studies  of  Brittle  Behaviour  of  Ceramic  Materials.  ASD  TR  61-628,  Part  II,  April  1962,  p.  38-42. 
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4.  Local  Stresses 


Loads  on  bend  specimens  applied  through  knife  edges  or  smal  1 -diameter  rollers 
result  in  high  stresses  under  these  line  loads.  High  compressive  contact  stresses 
can  result  and  cause  local  crushing.  Localized  contact  can  cause  a  more  subtle  prob¬ 
lem,  which  is  referred  to  as  the  wedging  stress. 

Also,  shear  stress  near  the  locality  of  the  load  point  can  be  several  times 
higher  than  that  predicted  by  beam  theory. 

a.  Contact  stress 

Reference  5  gives  equations  for  determining  the  contact  pressure  between  a 
cylinder  (or  roller)  and  a  flat  surface  (see  Figure  4)  as  a  function  of  the  applied 

load,  modulus  of  each  material,  and  the  roller  radius.  If  it  can  be  assumed  that  the 

two  materials  are  identical  and  that  the  allowable  bearing  pressure  or  contact  pres¬ 
sure  can  be  as  high  as  twice  the  bend  strength  of  the  material,  then  limits  on  the 
roller  radius  for  both  loading  systems  will  result.  For  example,  from  Reference  5 
we  have: 

P(nax  =  0.59/PE/2bp^  (17) 

where  pmax  is  the  maximum  contact  pressure.  (Note  that  the  roller  radius  can  be 
either  px  or  p2.)  However,  we  shall  assume  that  p^x  —  2oj>.  Also  for  four-point 

loading,  05  =  6Pa/bd2,  and  for  three-point  loading,  =  (3/2)  PL/bd2 .  Substituting 
of  ob  into  (17)  and  solving  for  pi/d  we  obtain 

Pi/d  2  7.25/a/d,  for  the  four-point  loaded  beam,  and  (17a) 

Pi/d  2  29.0/L/d,  for  the  three-point  loaded  beam,  (17b) 

where  it  was  assumed  that  E/o^  =  1  x  103. 


P  p  Figure  4.  Contact  point  tangency  shift. 


b.  Wedging  stress 

The  effect  of  the  wedging  stress  is  to  provide  a  substantial  tensile  stress  con¬ 
tribution  at  the  compressive  side  of  the  beam  adjacent  to  the  load  points.  A  net  ten¬ 
sile  stress  can  not  be  created  if  d/2  V  <  1,  according  to  Reference  14.  More  impor¬ 
tantly,  a  tensile  stress  is  added  to  that  already  present  due  to  beam  bending  at  the 
tensile  side  of  the  beam,  thereby  causing  a  deviation  from  the  assumed  stress  calcu¬ 
lated  by  simple  beam  theory. 

This  problem  is  generally  treated  in  Reference  2  and  particular  results  from 
von  Kdrmdn  and  Seewald18  for  a  similar  situation  are  used  to  estimate  this  error. 

An  analysis  for  this  error  is  given  in  Appendix  D.  The  resulting  error  determination 
for  four-  and  three-point  loaded  beams  is  given  in  Table  11.  In  the  calculation  of 
the  errors,  which  are  a  function  of  a/d  or  L/d,  as  well  as  x'/d»  the  computed  05 
corresponds  to  the  failure  site  location  (x'/d). 

c.  Beam  overhang 

The  overhangs  of  the  beam  must  be  great  enough  so  that  the  local  stresses  at 
the  beam  support  points  are  not  amplified  due  to  beam-end  effects.  These  stresses 
are  damped  out  within  a  distance  equal  to  one  beam  depth.18  Thus,  by  allowing 

Lt  >  L  +  2d  (18) 

we  avoid  beam-end  effects. 

d.  Fracture  origins 

Since  the  contact  stress  is  damped  out  with  a  distance  d,  fracture  origins 
for  the  four-point  loaded  beam  should  be  located  within  the  inner  span  and  be  no 
closer  than  the  beam  thickness  to  either  of  the  inner  loading  points  so  that  test 
data  results  can  be  considered  valid. 

The  distance  from  the  closest  support  point  to  the  fracture  origin  of  a  three- 
point  loaded  beam  should  be  used  to  determine  the  moment  and  thus  the  failure  stress. 
Also,  if  the  failure  origin  is  closer  to  either  of  the  support  points  rather  than  the 
center  load  point,  the  result  should  be  considered  invalid. 


5.  Contact  Point  Tangency  Shift 

Significant  changes  in  span  length  can  occur  in  both  four-point  and  three- 
point  loading  systems  if  contact  radii  of  support  and  load  points  are  large  com¬ 
pared  to  beam  depth.  The  shift  in  point  of  tangency,  as  shown  by  hx  and  h2  in 
Figure  4,  is  a  function  of  the  contact  radii,  specimen  thickness,  and  the  ratio  of 
the  modulus  of  elasticity  to  the  bend  strength.  For  materials  that  behave  elastic- 
ally,  such  as  those  considered  here,  the  change  in  tangency  point  and  thus  the  error 
arising  because  of  the  change  in  moment  arm  from  the  ideal  can  be  predicted  mathemati¬ 
cally  for  linear  systems.  This  is  accomplished  and  is  presented  in  Appendix  E.  The 


18.  VON  KARMAN.  T..  and  SEEWALD,  F.  Alhandi  Aerodynam,  Inat.  Tech.  Hochachule,  1946,  p.  256 
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approach  was  patterned  after  Westwater19  who  corrected  for  span  shortening  but  ignored 
friction  at  the  support  points  of  a  three-point  loaded  beam.* 

In  Appendix  E  the  formulas  are  derived  for  a  four-point  loaded  beam  and  then 
reduced  to  the  special  case  of  a  three-point  loaded  beam.  These  results  are  put  in 
terms  of  error  functions  assuming  the  simple  beam  theory  is  applied  without  correcting 
for  span  shortening,  as  in  the  case  of  the  lower  support,  and  span  lengthening  between 
the  upper  loading  points  shown  in  Figure  4. 

The  errors  are  determined  for  four-point  loaded  beams  of  1/3  and  1/4  loading 
points  as  a  function  of  pi/d  and  p2/d,  and  the  the  three-point  loaded  beam  as  a  func¬ 
tion  of  pi/d  only.  It  was  assumed  that  E/a=l*103.  These  errors  are  given  in  Table  12. 

6.  Surface  Preparation 


The  flexure  strength  of  each  brittle  material  is  not  only  supersensitive  to  the 
final  surface  finish  because  the  maximum  tensile  stress  occurs  at  the  beam  surface, 
but  is  also  highly  sensitive  to  prior  finish  history.  For  this  reason  it  is  impossi¬ 
ble  to  specify  an  optimum  surface  finish  procedure  for  all  brittle  materials,  so  that 
failure  will  be  due  to  inherent  flaws  related  to  the  material  or  material  processing, 
rather  than  an  imposed  defect  resulting  from  the  finish  process.  Indeed,  the  designer 
or  materials  developer  may  not  be  able  to  specify  a  particular  finish  procedure. 
Therefore,  rather  than  attempt  to  dictate  surface  finish  requirements,  it  is  suggested 
that  each  set  of  reported  test  data  results  be  accompanied  by  surface  finish  history 
and/or  material  process  history,  whichever  is  applicable. 


There  are,  however,  several  specific  recommendations  related  to  surface  finish¬ 
ing  procedures  that  can  be  presented.  Corner  flaws  resulting  from  chipping  or  crack¬ 
ing  during  the  grinding  operation  are  sources  of  low-strength  failure.  Rounding  or 
beveling  of  the  corner  as  depicted  in  Figure  5  appears  to  reduce  premature  failure.20 
Since  a  chamfer  will  double  the  number  of  edges,  thus  doubling  the  source  of  flaw 
locations,  rounding  is  preferred.21  Also,  it  is  important  to  grind  the  edges  and  flat 
surfaces21  by  a  motion  parallel  to,  rather  than  perpendicular  to,  the  specimen  length. 
It  is  further  indicated20  that  finishing  of  the  corner  should  be  comparable  in  all 
aspects  to  that  applied  to  the  beam  surfaces. 


Rounded  Corners  Chamfered  Corners 


Figure  5.  Beam  cross  section. 


‘Westwater  also  determined  an  approximate  relationship  for  the  horizontal  load  arising  because  of  tangency  shift.  However,  for 
beams  of  small  deflection,  the  error  is  negligible. 

19.  WESTWATER,  I.  W.  Flexure  Testing  of  Plastic  Materials.  Proc.  ASTM,  v.  49,  1949. 

20.  RICE,  R.  W.  Machining  of  Ceramics.  Proc.  of  the  Second  Army  Materials  Technology  Conference  -  Ceramics  for  High  Performance 
Applications,  J.  J.  Burke,  A.  E.  Gorum,  and  R.  N.  Katz,  ed..  Brook  Hill  Publishing  Company.  Chestnut  Hill,  Massachusetts,  1974. 

21.  RICE,  R.  W.  The  Effect  of  Grinding  Direction  on  the  Strength  of  Ceramics.  The  Science  of  Ceramics  Machining  and  Surface 
Finishing,  S.  J.  Scheider  and  R.  W.  Rice,  ed.,  NBS  Special  Publication  348,  Washington,  DC,  G.  P.  O.  (SD  Cat.  No.  13.10:348), 
1972,  p.  365-376. 
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If  the  corner  radii  or  chamfer  is  small,  the  error  in  ignoring  the  change  in  mo¬ 
ment  of  inertia  will  be  negligible.  The  limiting  ratio  of  comer  radii  or  45°chamfer 
dimension  to  beam  depth  can  be  determined  from  the  error  analysis  due  to  neglecting 
the  change  in  moment  of  inertia  given  in  Appendix  F.  This  error  in  determining  flexure 
stress,  when  neglecting  comer  radii  or  45°  chamfer,  is  given  in  Table  13.  Limiting 
ratio  of  comer  radii  or  45°  chamfer  is  indicated  by  a  line  at  an  error  level  of 
+0.5%  or  less. 


REFERENCE  STANDARD  BEAMS 


Geometry  Ratios 

From  previous  discussions  included  under  "Simple  Beam  Theory  Assumptions"  and 
"External  Influences"  a  range  of  practical  geometry  ratios  of  b/d,  a/d,  L/d,  and  pc/d 
can  be  chosen  so  that  the  error  span  is  tolerable.  This  is  given  in  Table  14,  as  well 
as  the  source  of  error.  Also  given  in  the  table  are  the  limits  on  roller  radius  to 
beam  depth  ratios  so  that  bearing  stress  fracture  at  the  contact  points  is  precluded. 
The  range  of  practical  geometry  ratios  are  5.0  <.  a/d  <  12.5  for  the  four-point  loaded 
beam  and  20  <  L/d  .s.  25  for  the  three-point  loaded  beam,  with  b/d  <.  15.0  and  pc/d  >  100. 
These  ratios  were  chosen  to  limit  the  accumulated  error  span  to  a  reasonable  total  of 
less  than  approximately  +2%.  A  discussion  of  errors  considering  the  external  influ¬ 
ences  as  well  will  be  subsequently  presented. 

Beam  Dimensions 

Most  workers  in  the  field  utilize  either  a  1/3  or  a  1/4-point  loaded  beam. 

For  convenience  we  shall  choose  L/a  =  3,  because  this  will  allow  the  four-point  beam 
length  to  be  equal  to  that  of  the  three-point  loaded  beams. 

We  have  already  specified  that  b/d  -i  15  and  we  shall  arbitrarily  let  b/d  =  2.0 
and  L/d  =  21.0.  By  choosing  a  convenient  and  practical  beam  width,  such  as  0.250  in. 
(6.35  mm),  we  can  specify  all  other  beam  dimensions  for  both  types  of  reference  stan¬ 
dard  beams.  These  dimensions  are  shown  in  Table  15.  Notice  that  (7)  and  (8)  are  sat¬ 
isfied,  i.e.,  L/d  -  a/d  .s  150  and  L/d  i  150,  which  insures  negligible  nonlinear  effects 
due  to  friction  at  the  load  and  support  points. 

As  indicated  previously  it  is  not  the  intent  to  require  specific  dimensions  for 
the  two  beam  systems.  Indeed,  this  would  impose  a  considerable  waste  of  effort  be¬ 
cause  many  organizations  already  have  a  backlog  data  base  from  which  test  results  were 
previously  obtained  from  other  than  those  beam  dimensions  shown  in  Table  15.  Alterna¬ 
tively,  it  is  suggested  that  two  sets  of  beam  dimensions  indicated  in  the  table  be 
considered  as  reference  bases  to  which  other  beam  geometries  can  be  related  by  statis¬ 
tical  analyses.  The  most  generally  acceptable  analysis  appropriate  to  brittle  materi¬ 
als  is  that  conceived  by  Weibull.13  With  this  approach  in  mind,  strength  relation¬ 
ships  between  the  reference  standard  beams  and  beams  of  other  dimensions  for  volume- 
sensitive  and  surface-sensitive  materials  are  presented  in  the  following  section. 

Also  given  is  an  estimate  of  error  of  the  Weibull  parameters  as  a  function  of  sample 
size. 


STRENGTH  RELATIONSHIP  AND  SAMPLE  SIZE 

The  use  of  beam  configurations  other  than  those  recommended  in  Table  15  is  depen¬ 
dent  upon  the  establishment  of  a  failure-strength  relationship  between  the  nonstandard 
and  the  reference  standard  beam.  The  establishment  of  such  relationships  is  considered 
the  responsibility  of  each  investigator. 
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The  type  of  analysis  that  has  been  used  with  varying  degrees  of  success  to 
relate  failure  strengths  of  brittle  materials  is  attributed  to  Weibull.13  Many  in¬ 
vestigators  have  used  this  approach  to  relate  strength  levels  of  various  types  of 
specimen  configurations  either  on  a  stressed  volume  or  surface  area  basis.  5  The 
reader  is  cautioned  that  confirmation  of  such  an  analysis  or  lack  thereof  may  well 
depend  on  a  number  of  factors  including  the  test  material.  As  examples  of  such  cor¬ 
relation  and  lack  of  it,  Weibull  statistical  correlation  was  justified  by  Davis22 
for  reaction-bonded  silicon  nitride  but  inappropriate  for  Lewis'23  work  in  alumina 
fabricated  by  several  processes. 

A  computer  program  for  statistical  evaluation  of  composite  materials  is  available 
in  Reference  24.  This  program  determines  the  desirability  of  a  particular  probability 
density  function  in  predicting  fracture  strength  of  ranked  empirical  data.  The  can¬ 
didate  functions  include  normal,  log  normal,  and  Weibull.  Root  mean  square  error 
results  can  be  tabulated  for  each  functional  comparison.  The  effects  of  statistical 
ranking  can  be  readily  listed  in  the  computer  output. 

The  data  mean  and  standard  deviations  with  corresponding  levels  of  confidence 
can  be  included  in  the  printed  results.  The  Weibull  parameters,  obtained  from  the 
maximum  likelihood  method,  and  corresponding  confidence  intervals  can  be  obtained  from 
this  program. 

Since  a  Weibull-type  analysis  is  applicable  in  many  instances,  resulting  formulas 
for  the  simple  two-parameter  system25  to  determine  the  risk  of  rupture  for  the  four- 
point  and  three-point  loading  systems,  are  presented  below,  for  the  sake  of  complete¬ 
ness. 

Volume-Sensitive  Material 

If  the  strength  of  material  is  dependent  upon  the  volume  of  stressed  material, 
the  following  risk -of- rupture  equations  are  applicable  for  a  beam  stressed  in  pure 
bending*  (four-point  loading) : 

rrV  =  [V£/2(M+l)](ob/o0)»«  (19) 

and  for  a  beam  stressed  under  three-point  loading,  we  have: 

RRV  =  [VL/2(M+l)2](ob/aQ)M  ;  (20) 

where  R^y  is  the  risk  of  rupture  (volume  basis).  The  volume  for  a  four-point  loaded 
beam  in  (19)  is  =  tbd  and  the  volume  for  a  three-point  loaded  beam  in  (20)  is 
Vb  =  Lbd;  also,  ab  is  the  mean  fracture  stress;  a0  is  a  normalizing  parameter,  also 
called  the  characteristic  value;  and  M,  the  Weibull  slope,  also  called  the  shape  pa¬ 
rameter,  is  the  reciprocal  function  of  the  variability  of  the  failure  stresses  of  the 
specimens  in  the  sample.  Note  o0  has  units  of  stress  times  volume  (1/M). 

M  can  be  determined  by  a  number  of  different  methods  (see  References  22,  25,  26, 
and  27).  The  accuracy  by  which  M  can  be  determined  is  discussed  later  under  "Weibull 
Parameter  Estimate  and  Sample  Size." 

•Since  one  of  the  requisites  for  a  valid  test  is  that  failure  occur  within  the  center  span,  only  the  constant  moment  section  or  pure 
bending  form  is  considered  for  the  reference  standard. 

22.  DAVIS,  D.  G.  S.  The  Statistical  Approach  to  Engineering  Design  in  Ceramics.  Proc.  Br.  Ceramics  Soc.,  no.  22,  1973,  p.  429-452. 

23.  LEWIS,  D.,  and  OYLER,  S.  M.  An  Experimental  Test  of  Weibull  Scaling  Theory.  J.  Amer.  Cer.  Soc.,  v.  59,  no.  11-12,  November 
December  1976,  p.  507-510. 

24.  AMMRC  Interim  Report  No.  1  to  MIL-HDBK  17,  Composite  Materials  for  Aircraft  and  Aerospace  Applications  (Proposed),  December  1980. 

25.  DeSALVO,  G.  J.  Theory  and  Structural  Design  Applications  of  Weibull  Statistics.  Westinghouse  Astronuclear  Laboratory, 
WANL-TME-2688,  1970. 


If  the  average  strength  of  a  nonstandard  four-point  or  three-point  loaded  beam 
is  to  be  compared  to  that  of  its  reference  standard  beam,  the  following  relationships, 
obtained  by  simply  equating  the  risk  of  rupture,  are: 


gbNS4  _  /v'RS4\1/M  pbNS3  =  /VrS3\1/M 

abRS4  \VnS4/  ’  abRS3  \Vns3/ 


(21) 


where  pbNSi,,3  anc*  °bRS4  3  are  tbe  average  fracture  stress  in  bending  of  a  nonstandard 
and  reference  standard  fyeam,  and  3  and  Vrs4  are  the  volumes  of  the  four-point 

or  three-point  loaded  nonstandard  ana  reference  standard  beams,  i.e.,  Vrs4  =  £bd  or 


If  the  average  strength  of  the  reference  standard  three-point  loaded  beam  is  to 
be  compared  to  the  reference  standard  four-point  loaded  beam,  then 


pbRS3/abRS4  =  [(1/3)(M+1)]1/M; 


(22) 


note 


vRS4/vRS3  =  1/3- 


Surf  ace-  Sensi  tive  Material 

If  the  strength  of  the  material  is  dependent  upon  the  surface-stressed  material, 
the  following  general  risk-of -rupture  equations25  are  applicable  for  a  rectangular 
beam  stressed  by  pure  bending  (four-point  loading): 


rRS 


2(M+1) 


and  a  beam  stressed  by  three-point  loading: 


RRS 


— ik - (_L*  xVfb 

2 (X+l)  (M+l)  \M+1  /\a° 


(23) 


(24) 


where  Rpg  is  the  risk  of  rupture  (surface  basis);  X  is  defined  as  b/d;  in  (23)  is 
defined  as  2£(b+d)  ;  in  (24)  is  defined  as  2L(b+d) ;  and  Ms  is  the  Weibull  slope 
parameter  for  a  surface-sensitive  material. 


The  Weibull  slope  parameter  for  the  surface  material  ir.  general  is  not  the  same 
as  that  for  the  interior  material.  According  to  Paluszny  and  Wu27  these  differences 
can  be  minimized  by  improving  surface  finishes.  Regardless,  for  such  a  case,  (19) 
and  (20)  are  still  applicable. 


26.  McLFAN,  A,  F.,  and  FISHFR,  E.  A.  Brittle  Materials  Design.  High  Temperature  Gas  Turbine.  Ford  Motor  Company, 
Contract  DAAG46-7I-C-0I62,  Interim  Report,  AMMRC  CTR  77-20,  August  1977. 

27.  PALUSZNY,  A.,  and  WU,  W.  Probabilistic  Aspects  of  Designing  with  Ceramics.  Presented  at  the  22nd  Annual  Gas  Turbine 
Conference  of  A.S.M.F.,  Philadelphia.  Pennsylvania,  March  27-31,  1977. 
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However,  if  the  average  strength  of  a  nonstandard  four-point  or  three-point 
loaded  beam,  sensitive  to  surface  area,  is  to  be  compared  to  that  of  its  reference 
standard  beam,  the  following  relationships  are  appropriate: 

i^Vl/(M^1)>>«sdjI/Ms|  and 

%RS4  i\ANS4/\l/(M+l)^Ns4/)  ^ 

gbNS3  _  j/ARS3\  AnSs+A  /l/CM^j^RSsV)1^5 
gbRS3  (\ANS3/\XRS3+V  \l/(M+l)+XNs3/j 

where  Ans4  3  and  Ars4  3  are  surface  areas  of  the  four-point  or  three -point  loaded 
nonstandard  and  reference  standard  beams,  i.e.,  Apg^  =  2i(b+d)  and  Ars3  =  2L(b+d) 
where  Xmc  and  XDC  are  the  width-to-height  ratio  of  the  nonstandard  and  reference 
standard  Balms.  RS“’3 

Weibull  Parameter  Estimate  and  Sample  Size 

Flexure  tests  on  hot-pressed  silicon  nitride  material  reported  by  McLean  and 
Baker28  show  the  effect  of  Weibull  slope  M  for  specific  component  reliability.  The 
strength  requirement  for  a  specific  component  reliability  was  decreased  16%  by  a 
reported  20%  increase  in  M  from  a  nominal  value  of  10,  and  was  increased  27%  by  a 
20%  decrease  in  the  slope. 

Different  techniques  will  produce  somewhat  different  results,  according  to 
McLean  and  Fisher,26  when  estimating  the  Weibull  parameters.  Two  statistical  methods 
had  been  used  during  preliminary  analysis  of  hot-pressed  silicon  nitride  material 
strength  data,  and  results  indicated  that  the  estimate  of  the  characteristic  value 
o0  (or  scale  parameter)  were  very  close  while  the  Weibull  slope  estimates  vary  and 
thus  would  yield  considerable  differences  in  the  component  strength  requirement. 

The  following  is  quoted  directly  from  Reference  26  (except  to  change  reference 
and  figure  numbers  appropriate  for  this  report)  because  it  succinctly  addresses  the 
answer  to  the  question  of  proper  sample  size:  "The  exact  confidence  intervals  for 
the  parameters  are  based  on  the  distributions  obtained  by  Monte  Carlo  methods  presented 
in  Thoman  et  al.29  It  is  not  unexpected  that  the  uncertainty  in  the  estimation  of  a 
parameter  will  increase  as  the  sample  size  decreases.  This  uncertainty,  however,  has 
rarely  been  quantified.  The  width  of  the  confidence  intervals  for  the  parameters  is 
a  measure  of  the  uncertainty  and  aids  in  the  selection  of  the  sample  size  of  a  test. 
Figures  6  and  7  are  drawn  from  Reference  29  and  show  the  90%  confidence  bounds  for 
the  Weibull  slope  and  the  characteristic  value."  (Figure  7  differs  from  that  given 
in  Reference  26  in  that  two  additional  M  values  were  computed  and  shown.)  "The  bounds 
for  the  Weibull  slope  are  a  function  of  sample  size  only,  while  for  the  characteristic 
value  they  are  a  function  of  both  the  sample  size  and  the  Weibull  slope.  As  can  be 
seen  from  the  graphs,  the  error  or  uncertainty  in  estimates  from  small  sample  sizes 
is  very  large.  Important  judgements  and  significant  analysis  should  not  be  based  on 
small  samples.  Sample  sizes  of  at  least  30  should  be  used  for  all  but  the  most  pre¬ 
liminary  investigations.  An  uncertainty  of  ±10%  in  Weibull  slope  requires  more  than 
120  samples.  This  uncertainty  is  not  peculiar  to  just  ceramics,  but  is  intrinsic  to 
the  statistical  analysis  of  data,  whether  that  data  be  material  strength  or  the  life 

28.  McLEAN,  A.  F.,  and  BAKER,  R.  R.  Brittle  Materials  Design,  High  Temperature  Gas  Turbine.  Ford  Motor  Company. 

Contract  DAAG46-71-C-0162,  Interim  Report,  AMMRC  CTR  76-31 ,  October  1976. 

29.  THOMAN,  D.  R.,  BAIN,  L.  J.,  and  ANTLE,  C.  E.  Inferences  on  the  Parameters  of  Weibull  Distribution.  Technometrics, 
v.  11,  August  1969,  p.  445-460. 
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Characteristic  Value  Error  Percent  of  Expected  Value 


of  some  electronic  component.  The  choice  of  sample  size  depends  on  many  factors  in¬ 
cluding  the  cost  and  timing  of  testing  and  the  degree  of  conservation  which  is  accept¬ 
able,  but  erroneous  judgements  may  be  made  and  unacceptable  designs  pursued  if  the 
sample  sizes  are  too  small.” 

LOADING  SPEED 

It  is  well  known  that  brittle  materials  are  strain-rate  sensitive,  and  thus 
speed  of  loading  will  influence  the  stress  at  which  failure  of  the  beam  will  occur. 

To  choose  a  "static”  speed  which  would  insure  no  strain-rate  effect  for  all  brittle 
materials  would  not  only  be  impossible  but  time  consuming  when  testing  those  materials 
that  exhibit  a  lesser  strain-rate  effect.  However,  most  materials  testing  facilities 
utilize  an  Instron  testing  machine.  Normally  the  slowest  test  speed  in  such  a  ma¬ 
chine  is  0.02  in. /min  (0.51  mm/min).  It  would  seem  appropriate,  therefore,  to  utilize 
this  speed  of  testing  in  conjunction  with  the  two  reference  standard  beams  given  in 
Table  15  to  establish  a  reference  standard  strain-rate.  This  is  accomplished  in  the 
following  for  both  the  1/3-four-point  and  three-point  loaded  beams. 

The  strain  rate  is  defined  as: 

e  =  (ob/E)/t  ,  (26) 

where  t  is  the  time  of  the  applied  load;  but  since  the  speed  of  loading  is  s=y/t,  then 

i  =  (°b/E)/ys  (27) 

where  y  is  the  deflection  of  the  beam  and  s  is  the  constant  speed  of  the  testing 
machine. 

The  deflection  at  the  center  of  a  four-point  loaded  beam  is: 

y  =  (Pa/24EI) (3L2-4a2)  .  (28) 

Recalling  that  Pa  =  2crbI/d,  and  substitution  of  this  and  (28)  into  (27)  gives: 

.  _  12ds  (29) 

3L2-4a2 

for  the  four-point  loaded  beam. 

For  the  1/3-four-point  loaded  beam,  where  a=L/3,  we  obtain: 

e  =  (108/23) (d/L2)s  .  (30) 

Using  the  same  approach  as  above,  we  obtain  the  strain  rate  for  the  three-point 
loaded  beam  as: 

i  =  (6d/L2)s  .  (31) 

With  a  loading  speed  of  0.02  in. /min  (0.51  mm/min)  or  333  x  10  6  in. /sec 
(8458  x  10’6  mm/sec) ,  we  obtain  strain  rates  of  28  x  10  6/sec  for  the  1/3-four-point 
reference  standard  beam  (29)  and  36  x  10~6/sec  for  the  three-point  reference  standard 
beam  (31).  These  strain  rates  shall  be  considered  as  reference  standards  for  both 
beam  systems  as  indicated  in  Table  15. 
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DISCUSSION  OF  ERRORS 


The  errors  associated  with  flexure  testing  of  beams  can  be  classified  according 
to  their  source,  i.e.,  related  to  the  use  of  simple  beam  theory  (beam  dimensions)  or 
sources  arising  from  external  load  applications.  Table  16  lists  the  sources  and  sum¬ 
marizes  the  error  span  for  the  first  category  and  Table  17  for  the  second. 

Table  16  shows  resulting  error  spans  of  -0.4%  to  +0.3%  for  both  the  1/3-four- 
point  reference  standard  beam  and  the  three-point  reference  standard  beam  geometries 
described  in  Table  15.  These  errors,  as  expected,  are  rather  small. 

The  errors  summed  in  Table  17  are  not  of  small  magnitude  and  require  some  ex¬ 
planation.  These  are  discussed  in  the  order  in  which  they  appear. 

Error  Due  to  Load  Mislocation 

In  the  calculation  of  the  error  associated  with  load  mislocation,  it  was  assumed 
that  loading  heads  and  supports  are  rigid,  can  not  rotate  (see  Figures  lb  and  2b), 
and  are  located  manually.  Thus  it  is  also  assumed  that  the  best  possible  location 
tolerance  that  can  be  obtained  is  ±1/64  in.  (0.40  mm).  The  errors  due  to  load  mis¬ 
location  were  determined  by  interpolation  from  Tables  6  and  8  with  e/L  =  ±0.006. 

Error  Due  to  Beam  Twisting 

Again  it  is  assumed  that  the  loading  heads  and  supports  are  designed  as  indica¬ 
ted  in  Figures  lb  and  2b.  Thus,  twisting  of  the  specimen  can  occur  due  to  nonparallel 
pairs  of  load  or  support  points  and  specimen  skewness,  as  schematically  indicated  in 
Figure  3.  Table  10  allows  the  error  determination  caused  by  such  undesirable  action. 
We  shall  determine  the  error  for  both  reference  standard  specimens  given  in  Table  16 
with  an  assumed  but  realistic  set  of  test  parameters  and  test  results  for  illustrative 
purposes . 

It  will  be  assumed  that  the  material  is  hot-pressed  silicon  nitride  with  E  = 

45  x  10G  psi,  and  v  =  0.25.  The  bend  strength  determined  during  the  test  for  the 
1/3-four-point  loaded  beam  was  60  x  103  psi,  and  for  the  three-point  loaded  beam  was 
67  x  103  psi.  The  measured  angle  of  twist  of  the  specimen  configurations  is  assumed 
to  be  small  with  respect  to  that  of  the  loading  fixtures.  It  is  assumed  that  the 
best  alignment  of  the  loading  fixtures  with  each  other  will  be  no  less  than  1.0°; 
thus  <j>s  s  0  and  <J>p  =  1°. 

In  order  to  utilize  Table  10,  we  must  first  determine  n;  recalling  equation  14b: 

n  =  [3ki  (E/ ab)/  (1+v)  ]  [  (d/LT)4>s  +  (d/  £ ’ ) $p]  (b/ £ ' )  . 

Since  b/d  =  2.0,  then  ki  from  Table  9  is  0.229;  for  the  1/3-four-point  loaded  beam, 
d/Lj  =  0.043,  d/£ '  =  0.143,  and  b/£ '  =  0.286.  Thus,  n  according  to  the  above  rela¬ 
tionship  is  approximately  0.4.  Utilizing  Table  10b,  with  b/d  =  2.0  and  £ ' /b  =  3.50, 
by  linear  interpolation  we  determine  the  error  to  be  -1.0%.  For  the  three-point 
loaded  beam,  d/Lj  =  0.043,  d/£'  =  0.095,  and  b/£ '  =  0.190;  and  n,  calculated  from 
the  above  relationship,  is  approximately  equal  to  0.2.  Entering  Table  10a  and  inter¬ 
polating  with  b/d  =  2.0  and  i'/b  =  5.25,  we  determine  the  error  to  be  approximately 
-0.1%. 


Note  that  in  the  above  example,  n  <  1.0,  indicating  that  failure  occurs  after 
bottoming.  If  n  had  been  greater  than  1.0,  i.e.,  no  bottoming  prior  to  failure,  then 
Table  lOe  would  have  been  used  to  calculate  the  error. 


Error  Due  to  Friction 


If  moveable  roller  contact  points  are  not  used  in  the  load  and  support  fixtures, 
an  error  in  determining  the  flexure  stress  at  failure  will  result.  An  estimate  of 
the  coefficient  of  friction  obtained16  by  comparing  fracture  strengths  when  using 
rigid  as  compared  to  roller-type  contacts  was  determined  to  be  0.4  for  silicon  ni¬ 
tride  material.  This  value  was  also  obtained  in  Reference  15  for  hard  steel  on  gar¬ 
net  and  Reference  17  for  steel  on  graphite.  Using  equations  15  and  16,  the  errors  as 
indicated  in  Table  17  are  +6.1%  for  the  four-point  beam  and  +4.0%  for  the  three-point 
beam. 

Error  Due  to  Wedging  Stress 

Since  a/d  for  the  four-point  reference  standard  beam  is  7  and  L/d  is  21  for  the 
three-point  reference  standard  beam,  then,  from  Tables  11a  and  lib,  assuming  the  max¬ 
imum  error  can  occur,  i.e.,  x'  =  0,  the  respective  errors  are  +0.7%  and  +0.9%. 

Error  Due  to  Contact  Point  Tangency  Shift 

There  are  two  conflicting  requirements  regarding  contact  radius  at  the  loading 
and  support  points:  the  first  is  that  radii  must  be  great  enough  so  that  contact  or 
bearing  pressure  does  not  cause  local  failure  of  the  beam,  and  the  second  is  that  the 
contact  radii  be  small  enough  so  that  the  error  due  to  contact  point  tangency  shift 
is  not  great. 

Equations  17a  and  17b  allow  the  determination  of  the  radii  such  that  the  contact 
pressure  is  not  excessive.  Since  we  have  allowed  the  contact  radius  to  beam  thickness 
ratio  to  be  1.5  (pi/d  =  p2/d  =  p/d  =  1.5)  for  both  reference  standard  beam  systems, 
equations  17a  and  17b  are  both  satisfied,  and  the  resulting  radius  of  3/16  inch 
(4.76  mm)  is  of  practical  size.  The  errors  due  to  contact  point  tangency  shift,  +0.5% 
and  +0.2%  were  then  obtained  from  Tables  12a  and  12c  for  the  four-point  and  three- 
point  beam  systems. 

Corner  Radii 

The  error  when  ignoring  corner  radii  in  determining  the  flexure  stress  at  fail¬ 
ure  and  neglecting  the  change  in  the  moment  of  inertia  with  b/d  =  2.0  and  r/d  =  0.06 

is  +0.4%  for  both  systems  as  obtained  from  Table  13. 

It  is  noted  that  three  error  sources  associated  with  the  determination  of  flex¬ 

ure  stress  at  failure,  not  included  in  Table  17  because  they  are  intrinsic  to  the 
material,  are  unequal  Young's  moduli  (Table  2),  anisotropy,  and  nonhomogeneity.  Each 
investigator  will  have  to  establish  if  these  factors  are  present,  and  if  so,  their 
magnitudes. 

Error  Span  Summary 

The  error  spans  indicated  in  Table  17  for  the  four-point  and  three-point  beam 
reference  standard  systems  are  -5.0%  to  +7.7%  and  -0.1%  to  +  5.6%.  Adding  those  given 
in  Table  16  results  in  a  total  of  -5.4%  to  +8.0%  and  -0.5%  to  +5.9%.  The  three-point 
loaded  beam  has  a  lower  error  span  associated  with  it  than  the  four-point  loaded  beam 
because  of  lesser  inherent  errors  in  load  mislocation,  beam  twisting,  friction  at  the 
load  and  support  contact  points,  and  contact  point  tangency  shift.  For  the  four-point 
loaded  beam,  the  two  major  sources  of  error  are  caused  by  load  mislocation  and  fric¬ 
tion,  which  can  be  reduced  to  a  great  degree  by  proper  design.  If  these  two  sources 
of  error  could  be  minimized  (assumed  to  be  zero),  the  total  error  spans  would  be 
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reduced  to  -1.0%  to  +1.6%.  The  major  source  of  error  for  the  three-point  beam  in 
bending  is  due  to  friction  and  if  eliminated  would  result  in  a  total  error  span  of 
-0.1%  to  +1.6%. 

These  minimized  total  error  spans  are  attainable  by  designing  load  and  support 
fixtures  that  incorporate  rolling  contacts  and  pivot  in  the  two  required  directions, 

i.e.,  about  the  beam  axis  and  transverse  to  it.  Such  loading  systems  have  been  de¬ 
signed14  and  extensively  used  in  practice.* 

RECOMMENDATIONS  AND  CONCLUSIONS 

1.  The  recommended  beam  configurations,  called  reference  standards,  are 
1/3-four-point  and  three-point  loaded  beams  having  the  following  dimensions  (see 
Figures  la  and  2a) : 

b  =  0.2S0  inch  (6.35  mm) 
d  =  0.125  inch  (3.18  mm) 
a  =  0.875  inch  (22.23  mm) 

L  =  2.625  inch  (66.68  mm) 

L-p  >  2.875  inch  (73.03  mm) 
r  >_  0.008  inch  (0.20  mm) 

2.  The  recommended  contact  radius  used  in  both  the  loading  and  support  fixtures 
is  0.1875  inch  (4.76  mm). 

3.  The  errors  due  to  simple  beam  theory  assumptions  for  the  reference  standard 
beams  are  relatively  small. 

4.  The  total  estimated  error  spans  for  the  above  beam  dimensions  are  -5.4%  to 
+8.0%  for  the  1/3-four-point  beam  and  -0.5%  to  +5.9%  for  the  three-point  beam.  Not 
included  in  these  error  spans  are  the  intrinsic  effects  of  unequal  Young's  modulus, 
anisotropy,  and  nonhomogeneity  of  the  beam  material. 

5.  If  a  properly  designed  fixture  is  utilized  during  testing  of  the  four-point 
loaded  reference  standard  beam,  such  that  errors  due  to  an  assumed  ±l/64-inch  load 
mislocation  and  an  assumed  1°  twisting  angle  of  the  loading  fixtures  are  eliminated, 
as  well  as  eliminating  a  coefficient  of  friction  of  0.4,  then  the  error  span  is  re¬ 
duced  accordingly  from  -1.0%  to  +1.6%. 

6.  The  major  source  of  error  associated  with  the  reference  standard  three- 
point  loaded  beam  is  due  to  friction  at  the  support  and  load  points;  elimination  of 
this  results  in  an  error  span  of  -0.1%  to  +1.6%. 

7.  It  is  recommended  that  the  load  and  supporting  fixtures  for  both  reference 
standard  systems  be  designed  so  that : 

a.  the  load  and  support  points  be  accurately  fixed  and  known;  and 

b.  friction  be  minimized  at  the  load  and  support  points  through  the 
use  of  freely  pivoting  rollers  or  by  other  means. 

•Private  communications  with  L.  Gibson  of  Carborundum  Co..  W.  Duckworth  of  Battelle  Columbus  Laboratories,  and  J.  Caverly  of 
Ford  Motor  Company. 


8.  In  addition  to  item  7  above,  the  reference  standard  four-point-beam  system 
employs  load  fixtures  that  are  free  to  pivot  in  two  mutually  perpendicular  directions. 
Although  two  beam  systems  are  recommended  and  referred  to  as  reference  standards, 
other  systems  are  acceptable  as  long  as  strength  relationships  through  proven  statis¬ 
tical  analyses  are  insured,  and  the  errors  associated  with  the  beam  and  loading  geom¬ 
etry  (obtainable  from  the  tables  in  the  text)  are  estimated  and  reported  along  with 
the  strength  data. 

9.  A  minimum  reference  standard  sample  size  of  30  is  recommended.  This  will 
result  in  an  error  of  approximately  ±25%  in  the  determination  of  the  Weibull  slope 
parameter  M  within  a  90%  confidence  band  (see  Figure  6).  Such  a  sample  size  will 
result  in  error,  dependent  on  M,  obtained  from  Figure  7  in  the  determination  of  the 
characteristic  value  cr0;  this  error  should  be  reported  also.  If  the  minimum  sample 
size  can  not  be  realized,  the  errors  in  determining  M  and  oQ  should  be  reported 
along  with  the  test  results. 

10.  No  specific  surface  finishing  procedure  is  recommended  but  this  will  be  con¬ 
sidered  the  responsibility  of  the  individual  investigator.  Nevertheless,  each  set 

of  test  results  should  include  surface  finish  or  material  process  history,  which¬ 
ever  is  applicable.  Surface  preparation  parallel  to  the  beam's  long  axis  is  recom¬ 
mended,  and  should  include  rounding  or  beveling  of  the  beam  edges.  Rounding  of  the 
edges  is  preferred,  but  both  are  acceptable. 

11.  If  the  fracture  origins  in  four-point  loaded  beams  can  be  determined,  they 
must  be  located  within  the  inner  span  and  be  no  closer  than  the  beam  thickness  to  the 
load  points  in  order  that  test  data  results  be  considered  valid. 

12.  The  distance  from  the  closest  support  point  to  the  fracture  origin  of  a 
three-point  loaded  beam  must  be  measured  and  used  in  computing  the  moment  and  result¬ 
ing  bending  stress.  Also,  if  the  failure  origin  is  closer  to  either  of  the  outer 
support  points  than  the  central  loading  point,  the  results  can  be  considered  invalid. 

13.  The  reference  standard  test  speed  of  loading  should  be  0.02  inch/minute 
(0.52  mm/min) .  This  results  in  a  strain  rate  ofg28  x  10  6/sec  for  the  reference 
standard  1/3-four-point  beam  system  and  36  x  10  /sec  for  the  three-point  beam  sys¬ 
tems.  If  any  other  speed  of  loading  is  used,  it  should  be  reported. 

14.  It  is  further  recommended  that  accurate  analyses  of  the  following  problems 
as  related  to  brittle  materials  be  accomplished: 

a.  Exact  stress  distribution  for  four-point  and  three-point  loaded  beams 
so  that  a  more  accurate  error  analysis  than  that  given  in  Table  1 

be  realized. 

b.  The  stress  distribution  for  parabolic  nonhomogeneous  material  for  four- 
point  and  three-point  loaded  beams  such  that  realistic  errors  may  be 
obtained  for  a  varying  modulus  of  elasticity  as  a  function  of  beam  depth. 

c.  Anisotropy  be  accounted  for  when  determining  bend  strength. 

15.  Finally,  it  is  suggested  that  a  simple  and  easy-to-use  loading  fixture,  as 
suggested  in  recommendations  7  and  8  above,  be  developed,  if  feasible. 


24 


TABULATIONS  OF  ERRORS  IN  CALCULATING  FLEXURE  STRESS 


Unless  otherwise  stated,  the  percent  error  is  determined  throughout  the  text  as 
[(°b-°x)/ax]100;  where  cr5=6M/bd2  and  ctx  is  more  nearly  the  true  bending  stress. 


Table  1.  ERROR  DUE  TO  NON-LINEAR 

Table  4. 

ERROR  CAUSED  BY  EFFECT 

STRESS  DISTRIBUTION 

OF  ANTICLASTIC  CURVATURE 

a/d  or  L/d 

%  Error 

E/ob  =  1*103 

0 

100 

b/d 

%  Error 

0.5 

51.6 

1.0 

0 

0.75 

32.2 

15.0 

0 

1.0 

21.1 

20.0 

0.1 

1.5 

10.6 

30.0 

0.6 

2.0 

6.2 

40.0 

1.5 

2.5 

4.1 

50.0 

2.6 

3.0 

2.9 

100.0 

4.7 

3.5 

2.1 

500.0 

5.9 

4.0 

1.6 

1000.0 

6.1 

4.5 

1.3 

co 

(-v2)100  =  -6.25% 

5.0 

1.1 

6.0 

0.7 

Note:  All 

errors  are  negative. 

8.0 

0.4 

10.0 

0.3 

00 

0 

Note:  All  errors  are  negative. 

Table  5.  ERROR  FOR  BEAMS  WITH  LARGE  DEFLECTION 

Table  2.  ERROR 

WHEN  Et  /  Ec 

£/  °b -  0.5  x  103 

ET/EC  %  Error 

Et/Ec  %  Error 

%  Error 

Beam  Loading 

0.20  +38.2 

1.025  -0.6 

0.40  +22.5 

1.050  -1.2 

L/d 

1/3-Four-Point  1/4-Four-Point  Three-Point 

0.60  +12.7 

1.075  -1.8 

0.80  +5.6 

1.10  -2.4 

0 

0 

0  0 

0.90  +2.6 

1.15  -3.5 

25 

0.4 

0.5  0.3 

0.925  +1-9 

1.20  -4.6 

50 

0.9 

1.1  0.7 

0.950  +1.3 

1.30  -6.5 

100 

2.0 

2.3  1.5 

0.975  +0.6 

1.40  -8.4 

150 

2.8 

3.2  2.1 

1.00  0 

1.60  -11.7 

250 

6.3 

7.1  4.7 

1.80  -14.6 

2.0  -17.2 

Note : 

All  errors 

are  positive. 

Table  3.  ERROR  CAUSED  BY 

INITIAL  BEAM  CURVATURE 

Table  6.  ERROR  DUE  TO  ECCENTRIC  LOAD  APPLICATION 
FOR  A  1/3-FOUR-POINT  LOADED  BEAM 

When  1/ L  =  1/3  and  aj/L  +  aj/L 

Vd 

%  Error 

e/L  =  ±(a,/L  -  1/3) 

%  Error 

1 

35.1 

0 

0 

2 

16.7 

0.0019 

1.0 

3 

10.9 

0.0038 

2.6 

4 

8.4 

0.0057 

3.8 

10 

3^2 

0.0076 

4.9 

15 

2.2 

0.0095 

6.0 

20 

1.7 

0.0114 

7.0 

40 

0.8 

0.0133 

8.1 

100 

0.3 

0.0333 

16.1 

0.0433 

18.7 

t  ~ 

100  [(ab"“c  )/“<;] 

0.0443 

18.9 

Note:  All 

errors  are  negative. 

Note:  All  errors  are  negative. 
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Table  7.  ERROR  DUE  TO  ECCENTRIC 
LOAD  APPLICATION  FOR  A 
1/4-FOUR-POINT  LOADED  BEAM 


Table  8.  ERROR  DUE  TO  ECCENTRIC  LOAD 
APPLICATION  FOR  A  THREE-POINT 
LOADED  BEAM 


When  t/ L  *  1/2  and  at/L  t  a2/L 

When  a^L  i  a2/L  /  1/2 
e/L  =  1/2  -  aj/L 

Table  9.  k) 

AND  k2 

e/L  =  t(aj/L  -  1/*) 

*  Error 

±e/L 

%  Error 

b/d 

kI 

^2 

0 

0 

0 

0 

1.0 

0.1406 

0.208 

0.0040 

3.8 

0.025 

0.25 

1.2 

0.166 

0.219 

0.0080 

7.1 

0.050 

1.0 

1.5 

0.196 

0.231 

0.0120 

10.0 

0.075 

2.3 

2.0 

0.229 

0.246 

0.0160 

12.6 

0.100 

4.2 

2.5 

0.249 

0.258 

0.0200 

14.7 

0.150 

9.9 

5.0 

0.291 

0.291 

0.0240 

16.6 

0.200 

19.0 

10.0 

0.312 

0.312 

0.0280 

18.3 

0.250 

33.3 

00 

0.333 

0.333 

0.0320 

19.8 

0.300 

56.3 

0.0340 

20.8 

0.400 

177.8 

0.0400 

21.8 

0.450 

426.3 

0.0465 

22.9 

0.500 

co 

Note:  All  errors  are  negative.  Note:  A11  errors  are  pos1tive- 


Table  10.  %  ERROR  DUE  TO  BEAM  TWISTING 

v  =  0.25 


I'/b 

b/d 

1.00 

1.2 

1.5 

2.0 

2.5 

5.0 

10.0 

« 

a.  n  *  0.20 

1.0 

3.18 

2.88 

2.61 

2.32 

2.12 

1.68 

1.47 

1.30 

2.0 

0.84 

0.76 

0.68 

0.60 

0.55 

0.43 

0.38 

0.33 

2.5 

0.54 

0.49 

0.44 

0.39 

0.35 

0.38 

0.24 

0.21 

5.0 

0.14 

0.12 

0.11 

0.10 

0.09 

0.07 

0.06 

0.05 

10.0 

0.03 

0.03 

0.03 

0.02 

0.02 

0.02 

0.02 

0.10 

00 

0 

0 

0 

0 

0 

0 

0 

0 

b.  n  =  0.40 

1.0 

10.58 

9.77 

8.94 

8.06 

7.44 

6.05 

5.36 

4.77 

2.0 

3.18 

2.89 

2.6! 

2.32 

2.12 

1.68 

1.47 

1.30 

2.5 

2.09 

1.89 

1.71 

1.51 

1.38 

1.09 

0.95 

0.84 

5.0 

0.54 

0.49 

0.44 

0.39 

0.35 

0.28 

0.24 

0.21 

10.0 

0.14 

0.12 

0.11 

0.10 

0.09 

0.07 

0.06 

0.05 

0 

0 

0 

0 

0 

0 

0 

0 

c.  n  =  0.60 

1.0 

19.01 

17.75 

16.51 

15.11 

14.11 

11  .79 

10.58 

9.54 

2.0 

6.58 

6.02 

5.48 

4. 10 

4.50 

3.61 

3.18 

2.81 

2.5 

4.44 

4.04 

3.67 

3.26 

2.99 

2.38 

2.09 

1.84 

5.0 

1.20 

1.08 

0.98 

0.86 

0.77 

0.62 

0.54 

0.48 

10.0 

0.31 

0.28 

0.25 

0.22 

0.20 

0.16 

0.14 

0.12 

oo 

0 

0 

0 

0 

0 

0 

0 

0 

d.  n  =  0.80 

1.0 

26.88 

25.38 

23.86 

22.12 

20.86 

17.92 

16.22 

14.79 

2.0 

10.58 

9.75 

8.94 

8.06 

7.44 

6.09 

5.36 

4.77 

2.5 

7.35 

6.73 

6.14 

5.50 

5.05 

4.09 

3.58 

3.17 

5.0 

2.09 

1 .89 

1.71 

1.51 

1.38 

1.10 

0.95 

0.84 

10.0 

0.54 

0.49 

0.44 

0.39 

0.35 

0.28 

0.24 

0.21 

« 

0 

0 

0 

0 

0 

0 

0 

0 

e .  n  =  1.00 

1.0 

33.75 

32.13 

30.47 

28.54 

27.11 

23.64 

21.73 

20.03 

2.0 

14.80 

13.74 

12.69 

11.53 

10.71 

8.83 

7.87 

7.05 

2.5 

10.58 

9.75 

8.94 

8.06 

7.44 

6.05 

5.36 

4.77 

5.0 

3.18 

2.89 

2.61 

2.32 

2.12 

1.68 

1.47 

1.30 

10.0 

0.84 

0.76 

0.68 

0.60 

0.55 

0.43 

0,38 

0.33 

00 

0 

0 

0 

0 

0 

0 

0 

0 

Note:  All  errors  are  negative. 
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Table  11.  ERROR  DUE  TO  WEDGING 


*'/d* 


Loading 

0 

0.125 

0.25 

0.375 

0.50 

0.75 

1.0 

1.50 

Four -point 

1.0 

4.7 

-0.5 

-2.8 

-2.1 

-1.4 

-0.7 

-0.3 

0 

1.5 

+3.1 

-0.3 

-1.9 

-1.4 

-0.9 

-0.5 

-0.2 

0 

2.0 

+2.3 

-0.2 

-1.4 

-1.1 

-0.7 

-0.4 

-0.2 

0 

3.0 

+1.5 

-0.2 

-1.0 

-0.7 

-0.5 

-0.2 

-0.1 

0 

4.0 

+1.1 

-0.1 

-0.7 

-0.5 

-0.3 

-0.2 

-0.1 

0 

5.0 

+0.9 

-0.1 

-0.6 

-0.4 

-0.3 

-0.1 

0 

0 

6.0 

+0.8 

-0.1 

-0.5 

-0.4 

-0.2 

-0.1 

0 

0 

8.0 

+0.6 

0 

-0.4 

-0.3 

-0.2 

-0.1 

0 

0 

10.0 

+0.4 

0 

-0.3 

-0.2 

-0.1 

-0.1 

0 

0 

15.0 

+0.3 

0 

-0.2 

-0.1 

-0.1 

0 

0 

0 

20.0 

+0.2 

0 

-0.1 

-0.1 

-0.1 

0 

0 

0 

40.0 

+0.1 

0 

-0.1 

-0.1 

0 

0 

0 

0 

60.0 

+0.1 

0 

0 

0 

0 

0 

0 

0 

CO 

0 

0 

0 

0 

0 

0 

0 

0 

Three-point 

L/d 

1.0 

+21.6 

-2.4 

-18.8 

-25.4 

-100.0 

+6.2 

+  1.3 

0 

1.5 

+13.4 

-1.4 

-10.4 

-10.2 

-10.1 

-100.0 

+2.6 

0 

2.0 

+9.7 

-1.0 

-7.2 

-6.4 

-5.3 

-5.5 

-100.0 

+0.1 

3.0 

+6.3 

-0.7 

-4.4 

-3.7 

-2.7 

-1.9 

-1.3 

-100.0 

4.0 

+4.7 

-0.5 

-3.2 

-2.6 

-1.8 

-1.2 

-0.6 

-0.1 

5.0 

+3.7 

-0.4 

-2.5 

-2.0 

-1.4 

-0.8 

-0.4 

0 

6.0 

+3.1 

-0.3 

-2.1 

-1.6 

-1.1 

-0.6 

-0.3 

0 

8.0 

+2.3 

-0.2 

-1.5 

-1.2 

-0.8 

-0.4 

-0.2 

0 

10.0 

+  1.8 

-0.2 

-1.2 

-0.9 

-0.6 

-0.3 

-0.2 

0 

15.0 

+  1.2 

-0.1 

-0.8 

-0.6 

-0.4 

-0.2 

-0.1 

0 

20.0 

+0.9 

-0.1 

-0.6 

-0.4 

-0.3 

-0.2 

-0.1 

0 

40.0 

+0.4 

0 

-0.3 

-0.2 

-0.1 

-0.1 

0 

0 

60.0 

+0.3 

0 

-0.2 

-0.1 

-0.1 

-0.1 

0 

0 

“ 

0 

0 

0 

0 

0 

0 

0 

0 

is  the  distance  on  either  side  of  the  load  contact  point  where  failure  occurs. 


Table  12.  %  ERROR  DUE  TO  TANGENCY  POINT  SHIFT 

E/ojj  =  1  x  101 


p2/d 


Loading 

pi/d 

1.0 

2.0 

5.0 

10.0 

a.  Four-point, 

1.0 

0.3 

0.4 

0.7 

1.2 

a/L  =  1/3 

2.0 

0.5 

0.6 

0.9 

1.4 

4.1 

0.9 

1.0 

1.3 

1.8 

6.1 

1 .3 

1.4 

1.7 

2.3 

8.2 

1.7 

1.8 

2.1 

2.7 

10.3 

2.1 

2.2 

2.6 

3.1 

b.  Four-point, 

0.67 

0.4 

0.6 

1.2 

2.2 

a/L  =  1/4 

1.35 

0.6 

0.8 

1.4 

2.5 

2.7 

1.0 

1.2 

1.8 

2.9 

4.1 

1.4 

1.6 

2.2 

3.3 

5.5 

1.8 

2.0 

2.6 

3.7 

6.9 

2.2 

2.4 

3.1 

4.1 

c.  Three-point, 

1.0 

0.1 

a/L  =  1/2 

2.0 

0.2 

4.0 

6.0 

0.4 
0.6  | 

Regardless  of 

n?/d 

value 

8.0 

0.8 

10.0 

1.0 

Note:  All  errors  are  positive. 
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Table  13.  %  ERROR  IN  DETERMINING  FLEXURE  STRESS 


a.  When  neglecting 

b/d 

corner  radii 

— ■ - 

r/d 

1.0 

2.0 

4.0 

0 

0 

0 

0 

0.02 

0.1 

0.1 

0 

0.04 

0.4 

0.2 

0.1 

0.06 

0.9  L 

0.4 

0.2 

0.08 

1.5 

0.8  I  0.4 

0.10 

2.4 

1.2 

0.6 

0.15 

5.4 

2.6 

1.3 

0.20 

9.7 

4.6 

2.2 

b.  When  neglecting 

c/d 

1.0 

2.0 

4.0 

45°  chamfer 

0 

0 

0 

0 

0.01 

0.1 

0.1 

0.1 

0.02 

0.2 

0.1 

0.1 

0.03 

JL5- 

0.3 

0.1 

0.04 

0.9  L 

0.5 

0.2 

0.05 

1 .4 

0.7  I 

I  0.4 

0.06 

2.0 

1.0  | 

!  0.5 

0.08 

3.4 

1.7 

0.9 

0.10 

5.2 

2.6 

1.3 

Note:  All  errors  are  positive. 


Table  14.  SOURCES  AND  ERRORS 


Four-Point 

Three 

-Point 

Error  Source 

a/d 

L 

p  i/d  r 

L/d  7 

pi/d  F 

Nonlinear  Stress 
Distribution,  (Table  1) 

>3.5 

i-2. 1 

- 

>7.0  <-0.6 

- 

Et  /  Ec,  (Table  2) 
Nonhomogeneity  and 
Anisotropy 

Error  Dependent  Upon  Material 

and/or  Fabrication  Process 

Initial  Beam 

Curvature,  (Table  3) 

pc/d  i  100, 

l  <  -0.3 

Anticlastic 

Curvature,*  (Table  4) 

b/d  -  15, 

=  0 

Large  Deflection, 

(Table  5) 

<12.5 

*0.3 

- 

i.25  <+0.3 

Wedging  Stress,** 

(Table  11) 

>5 

(Table 

<-0.9 

11a) 

- 

>20  i+O.9 

(Table  lib) 

- 

Tangency  Point 

Shift 

- 

- 

<1.5+  <+0.5 

(Table  12a,  a/L  =  1/3) 

<1.5  <+0.2 

(Table  12c) 

Neglecting  Corner 

Radii  or  Chamfer  (Table  13) 

r/d  <  0.06 
c/d  <  0.04 

t  <  +0.5 

Contact  Stress 

- 

- 

>.7.25/ (a/d) 

(Eg.  17a) 

>29.0/(L/d) 

(Eg.  17b) 

•Assumed  E/05  “  1*103 

••Assumed  the  most  conservative  case  (x ’ =  0). 
^Assumed  oJA  *  oj d  for  both  beam  systems. 
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Table  15.  REFERENCE  STANDARD  BEAM  SYSTEMS 


1/3-Four-Point  Beam 

Three-Point  Beam 

b/d 

2.0 

2.0 

a/d 

7.0 

10.5 

L/d 

21.0 

21.0 

vd 

1.5 

1.5 

d/Lr 

0.043  in. 

(1.104  mm) 

0.043  in. 

b 

0.250  in. 

(6.350  mm) 

0.250  in. 

d 

0.125  in. 

(3.175  mm) 

0.125  in. 

a 

0.875  in. 

(22.225  mm) 

1 .313  in.  (32.944  mm) 

r 

0.008  in. 

(0.203  mm) 

0.008  in. 

or  c 

0.005  in. 

(0.127  mm) 

0.005  in. 

i 

0.875  in. 

(22.225  ran) 

- 

L 

2.625  in. 

(66.675  mm) 

2.625  in. 

lt 

2.875  in. 

(73.025  mm) 

2.875  in. 

Pi 

0.1875  in. 

(4.763  mm) 

0.1875  in. 

s 

0.02  in. /min 

(0.508  mm/min)  0.02  in. /min 

e 

28x1 0~ 6 /sec 

36*10"6/sec 

Table  17.  %  ERRORS  CAUSED  BY  EXTERNAL  FACTORS 

FOR  REFERENCE  STANDARD  BEAMS 


Type  of  Loading 

Source 

1/3-Four-Point 

Three-Point 

Load  Mislocatlon 
e/L  =  ±0.006 

-4.0 

(Table  6) 

+0.1 

(Table  8) 

Beam  Twisting, 

Table  10 

(4s  =  4fr=  i°) 

-1.0  ' 

-0.1 

Friction 
(u  *  0.4) 

+6.1 

(Eq.  15) 

+4.0 

(Eq.  16) 

Wedging* 

+0.7 

+0.9 

(Table  11a) 

(Table  lib) 

Contact  Point 
Tangency  Shift, 

=  +0.5 

n  j/d  =  =  1.5 

(Table  12a) 

=  +0.2 
n,-  1.5 
(Table  12b) 

Corner  Radii, 
Table  13 

+0.4 

+0.4 

Error  Span 

-5.0%  to  +7.7'; 

-0. U  to  +5.6  . 

Assumed  the  most 

conservative  case,  x' 

=  0. 

Table  16.  %  ERRORS  DUE  TO  SIMPLE  BEAM  THEORY 

ASSUMPTIONS  APPLIED  TO 
REFERENCE  STANDARD  BEAM  TYPE  LOADING 


Source 

1/3-Four-Point  and 
Three-Point  Loading 

Nonlinear  Stress 
Distribution,  L/d=2 1 , 
Table  1 

-0.1 

Beam  Curvature, 
pc  >  100,  Table  3 

-0.3 

Anticlastic  Curvature, 
b/d  =  2,  Table  4 

0 

Large  Deflection,* 

L/d  =  21,  Table  5 

+0.3 

Error  Span  due  to 

Beam  Geometry 

-0.4  to  +0.3 

‘Assumed  E/c>b  *  0 . 5“ 10 3 
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APPENDIX  A.  ANTICLASTIC  CURVATURE 


When  a  beam  is  bent  by  a  moment,  it  produces  a  curvature  p  along  its  longitu¬ 
dinal  axis;  there  is  also  curvature  present  in  the  transverse  or  lateral  direc¬ 
tion.  This  moment  is  defined  by  orthodox  theory  as 

Mg  =  (EI/p)8  (A-l) 

where  6  is  a  parameter  representing  anticlastic  curvature  after  Ashwall.7 

Since 

ax  =  Mgc/I  =  (EI/pI)yS 

then 

ax  =  (Ey/p)g  (A-2) 

where  y  is  the  distance  from  the  neutral  axis. 


The  B  for  such  a  beam  of  the  simple  case  will  equal  1.0  and  if  the  beam 
width  to  depth  is  great,  i.e.,  b/d*-*”,  the  beam  can  be  considered  as  a  plate  so 
that  f?->-l/  (1-v2) .  It  is  worthwhile  to  know  the  intermediate  values  of  B  such  that 
the  effect  of  anticlastic  curvature  on  the  error  can  be  ascertained  when  assuming 
simple  beam  theory  (B  =  1.0)  is  valid. 


Ashwall7  has  determined  8  as  a  function  of  Poisson's  ratio,  beam  width, 
depth,  and  neutral  axis  curvature  by  accounting  for  anticlastic  curvature  and 
treating  the  structure  as  a  beam  on  an  elastic  foundation.  The  function  B  and 
related  terms  are  repeated  here  in  the  following: 


B  = 


1-v^ 


2  y8 


f(Yb) 


2/3 


yb/l-v2 


F  (yb) 


(A-3) 


where 


Y  = 


3(l-v2) 


d2p2 


f  (yb)  =  2(B2+C2) [sinh(yb)  +  sin(yb)] 
+  (B2-C2+2BC)  cosh(yb)  sin(yb) 
+  (B2-C2-2BC)  sinh(yb)  cos(yb) 
♦  2(B2-C2) (yb) , 


F(yb)  =  (B+C)  sinh(yb/2)  cos(yb/2) 

-  (B-C)  cosh(yb/2)  sin(yb/2) , 


*To  be  more  exact,  b2/pd  -*■  °°. 
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B  =  v  /sinh(yb/2)  cos(yb/2)  -  cosh(yb/2)  sin(yb/2)\  and 

"[TV)  \ 


/3(] 


C  = 


_v _  /si 

TV)  \ 


sinh(yb)  +  sin(yb)  ( 

inh(yb/2)  cos(yb/2)  +  cosh(yb/2)  sin(yb/2) 
sinh(yb)  +  sin(yb)  / 


The  calculations  performed  for  Table  4  in  the  text  were  accomplished  in  the 
following  manner: 

Since 


yb  =  b  a  /3 W  and  p  =  (E/o)v3,  then  substituting  into  the  above, 

V  d2p2 


allowing  y  =  d/2,  E/a  =  E/ajj  =  1  *  103, 


yb 


57.915  x  IQ"3  fb/d) 
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and  v  •-  0.25  for  ceramic  materials,  we  have; 


(A-4) 


By  programming  (A-3)  and  prescribing  b/d,  but  first  allowing  6=1.0,  then 
iterating  in  the  computer  through  (A-4),  the  relationship  between  b/d  and  6  was 
obtained.  Once  this  relationship  is  known,  the  percent  error,  defined  as 
l (l-6)/6] 100,  as  a  function  of  b/d  is  realized.  These  errors  are  given  in 
Table  4  as  a  function  of  b/d  with  E/a^  =  1  *  103  and  v  =  0.25. 


APPENDIX  B.  LOAD  MISL0CATI0N 


Four-Point  Loading 


Consider  the  usual  flexure  testing  setup  where  the  loading  head  is  rigidly 
attached  to  a  testing  machine,  schematically  shown  in  Figure  lb  and  idealized  in 
Figure  lc.  The  upper  loading  head,  where  the  inner  span  t  is  fixed,  can  only 
translate  in  the  vertical  direction,  and  the  lower  support  fixture,  where  the 
outer  span  L  is  fixed,  can  be  located  with  reference  to  the  loading  head.  The 
slope  and  deflection  equations  between  points  AB,  BC,  and  CD  are  as  follows; 


El  (dyAB/dx)  =  (Pxx2/2)  +  cxj 
El  yAB  =  (P i x3/ 6)  +  cxx  +  c2j 
El  (dygc/dx)  =  (P!-P2) (x2/2)  +  P2axx  +  c3 
El  yec  =  (P1-P2) (x3/6)  +  (P2aiX2/2)  +  c3x  +  c4 


(B-la) 

(B-lb) 


El  (dyCD/dx)  =  (Px-P2-P3) (x2/2)  +  P2axx  +  P3  (a x+f)  x  +  c5  ) 

[  (B-lc) 

El  yCD  =  (Pl-P2-P3)(x3/6)  -  (P2alX2/2)  +  P3  (ai<>  (x2/2)  -  c5x  +  c6) 
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where  Pi,  P2,  P3,  and  P4  are  the  loads,  x  and  y  are  defined  as  shown  in  Figure  1, 
E  is  the  Young's  modulus  of  the  material,  and  I  is  the  moment  of  inertia  of  the 
cross  section  of  the  beam. 

Through  the  use  of  the  various  boundary  conditions  the  constants  are  deter¬ 
mined  to  be: 


Ci  =  -1/6  {Pi (L2-a2)  +  (P2/L)  [(aj-L)3  +  (£+a2)a22]} 

C2  =  0 

C3  =  -1/6  (Pi(L2-a22)  +  (P2/L)  [3ai2L  -  (a 2+£)3  +  (£+a2)a22]} 

C4  =  P2aj3/6 

C5  =  -1/6  {P1L2  +  (P2/L)  [ai3+3aiL2-L3]  +  (P3/L)  [(aj+£)3  +  3(ai +£)  L2-L3]} 
C6  =  P3[(ai+£)3/6]  +  (P2ai3/6). 

In  order  to  determine  the  distribution  of  loading  between  the  vertical  loads 
Pi,  P2,  P3,  and  P4,  the  final  condition  of  equal  deflection  must  be  enforced  at 
locations  B  and  C  (see  Figure  lb) ,  which  is 


(yBc)x=aj  ^yBC-*x=ai+£. 

Enforcing  the  above  condition  in  the  second  part  of  Equation  B-lb  results  in 


Pl/P2  =  - (£/L) 2 


l/L  -  (l-ai/L)(2-^/L-2ai/L) 

(ai/L)3  -  (ai/L+£/L)3  +  {1  -  [l-(£/L)-ai/L]2}(i/L)  . 


(B-2) 


Utilizing  force  and  moment  equilibrium,  a  further  relationship  between  all  four 
forces  is  obtained  and  given  in  the  following: 

(Pl/P2) (ai/L)  +  (Pi/P2-1)£/L  P1/P2  +  ai/L  ‘  1 

P4/P3  =  - — — - - - - -  and  P3/P2  = -  (B-3) 

(Pl/P2)  -  1  +  ai/L  1  .  ai/L  _  £/L 

The  ratio  of  the  stress  at  x  (ox)  to  the  bending  stress  (ob)  from  Reference  14 
or  Equation  9  in  the  text,  where  it  is  assumed  that  ai  /  a2,  is: 


ax/ab  = 


P2+P3 


(B-4) 


where  a  is  the  value  at  ai  with  perfect  load  location  and  xi  is  defined 
as  shown  in  Figure  lc. 


By  manipulation  of  (B-2)  and  (B-3),  the  factor  Pi/(P2+  P3)  (B-4)  can  be  put 

into  terms  of  ai/L  and  l/L.  This  was  accomplished  and  the  results  are  shown  in  the 
following  equation: 
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0x'°b  = 


[ (ai/L) / (a/L) ] (1  -  ai/L  -  A/L) [ (A/L) (2  -  ai/L)-2(l  -  aj/L)2] 
3(ai/L)(l  -  A/L  -  ai/L)  -  (1  -  A/L)2 


By  defining  the  eccentricity  of  loading  as  e/L  =  ai/L  -  a/L,  Equation  B-5  becomes: 


ox/°b  = 


[(e/L  +  a/L) / (a/L) ] [1- (e/L  +  a/L)-  A/L] {  (A/L) [2- (e/L  +  a/L) ] - 2[1- (e/L  +  a/L)]2} 
3 (e/L  +  a/L)  [1  -  A/L  -  (e/L  +  a/L)]-(l-  A/L)2 


Calculations  of  ox/og  were  obtained  for  the  reference  standard  beams,  where  A/L  =  1/3 
as  well  as  A/L  =  1/2,  by  allowing  e/L  to  take  on  negative  values  only  in  Equation  B-6. 

Only  negative  values  were  considered  because  beam  failure  will  occur  due  to  a  realisti¬ 
cally  larger  moment  than  idealized  when  ignoring  eccentricity.  These  error  calculations 
although  determined  by  allowing  e/L  <  0  in  Equation  B-6,  are  indicated  as  ±e/L  in  Table  6 
for  A/L  =  1/3  and  Table  7  for  A/L  =  1/2.  This  simply  indicates  that  the  location  of  the 
maximum  moment  or  stress  is  at  Xi=  a2  when  e/L  <0  and  at  xj.=  ai  +  A  when  e/L  >  0. 


The  reader  is  cautioned  that  for  each  value  of  A/L  there  exists  a  set  of  limits  on 
(B-2),  (B-3),  and  (B-S).  That  is,  ax  can  be  such  that  either  P2  or  P3  can  equal  zero, 
because  (Ygc)x=ai  /  (^BC^x=ai+A-  anc*  system  changes  from  four-point  to  an  eccentric 
three-point  loading.  The  limiting  values  can  be  determined  by  allowing  P2  =  0  in  Equa¬ 
tion  B-2.  This  was  accomplished  in  terms  of  the  load  eccentricity  ratio  e/L  for  the  ref¬ 
erence  standard  beams,  when  A/L  =  1/3,  and  found  to  be  -0.0443,  and  when  A/L  =  1/2,  the 
result  was  -0.0465.  Thus  Tables  6  and  7  do  not  extend  beyond  e/L  of  -0.0443  and  -0.0465. 


APPENDIX  C.  BEAM  TWISTING 

if  line  loads  are  nonuniform  or  nonparallel  between  pairs  of  load  contacts, 
or  if  the  cross  section  of  the  specimen  is  skewed  along  its  length,  as  shown  in 
Figure  3,  a  net  torque  will  result.  The  addition  of  torque  gives  rise  to  a  max¬ 
imum  principal  stress  due  to  bending  and  torsional  stresses.  1,14  Failure  assumed 
to  be  caused  only  by  bending  stress  will  yield  an  error.  Two  cases  are  consi¬ 
dered:  Case  I  -  Failure  occurs  prior  to  specimen  realignment  in  the  bend  fixture 
(bottoming),  and  Case  II  -  Failure  occurs  at  or  after  bottoming. 

Case  I 

Recalling  that  the  bending  stress  for  a  loaded  beam  is 

ctx  =  ob  =  6Mfo/bd2  (C-l) 

where  Mg  is  the  measured  bending  moment  at  failure;  for  a  four-point  loaded  beam 
Mg  =  Pa,  and  for  a  three-point  loaded  beam  Mg  =  PL/4. 

The  shear  stress  due  to  torsion  is2 

t  =  Tb/k2bd2  (C-2) 

where  Tb  is  the  torque  and  equal  to  Pgb  for  four-point  bending  and  Pgb/2  for 
three-point  bending,  and  k2  is  a  numerical  factor  obtained  from  Reference  2  and 
is  given  in  Table  9. 
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Prior  to  bottoming,  the  normal  stress  is 

an  =  (crz/2)  (l-cos20)  +  (ax/2) (l+cos20)  +  rsin20. 


(C-3) 


where  0  is  the  angle  of  a  plane  inclined  to  the  x  axis.  Since  we  shall  assume  a 
plane  strain  condition,  i.e.  ,  ez  =  0,  then 

ez  =  0  =  (1/E) (az-vax)  or  az  =  vox  =  vcrb 

an  =  (ab/2)[(l+v)  +  (l-v)cos20]  +  isin20.  (C-4) 

Now  on  is  maximum  when 
tan20*  =  2x/(ax-oz); 
but  since 

o2  =  vax 

then 

tan20*  =  2t/(1-v)oij.  (C-5) 

where  0*  is  the  angle  of  a  plane  inclined  to  the  axis  at  which  principal  stress  is 
a  maximum.  Substitution  of  (C-l)  and  (C-2)  into  (C-5)  for  the  condition  prior  to 
bottoming  gives: 

tan20*  =  (Tb/Mb)/[3(l-v)k2],  (C-6) 

sin20*  =  (Tb/Mb)/[(Tb/Mb)2  +  (3k2 ) 2  Cl -v)  2 1*5,  (C-7) 

and 

cos20*  =  [3k2(l-v)]/[(Tb/Mb)2  +  (3k2)2(l-v)2]!S.  (C-8) 

Noting  that 

t  =  (ob/2)[(Tb/Mb)/3k2],  (C-9) 

and  by  substitution  of  the  above  relationships  into  (C-4)  with  some  algebraic  ma¬ 
nipulation,  we  obtain: 

°nmax  =  (<V2)  {fl+v)  +  (1/3k2)  [CV^)2  +  9(l-v) 2k22]*i }  (C-10) 

prior  to  bottoming.  The  shear  stress  due  to  torsion  can  be  related  to  the  twist 
angle2  of  the  beam  through  the  following  relationship: 

x  =  (k1/k2)Gd[(4>s/br)  +  (*F/*']  (C-ll) 

where  $s  is  the  twist  angle  along  the  length  of  the  specimen,  is  the  twist 
angle  between  a  pair  of  load  and  contact  points  relative  to  $s,  is  equal  to 
either  "a"  for  four-point  beam  systems  or  L/2  for  three-point  beam  systems,  kj 
is  another  numerical  factor2  given  in  Table  9  and  G  =  E/2(l+v),  the  shear  mod¬ 
ulus  of  the  material. 
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Equation  (C-2)  can  be  equated  to  (C-ll)  and  thus  we  obtain: 

Tbe  *  [kiE/2(l+v) ]bd2 [(d/LT)<j>s  +  (d/a')*p]  (C-12) 

where  Tb  j.s  the  torque  when  bottoming  occurs.  Thus,  in  order  for  (C-10)  to  be  ap- 
plicablef  Tb  must  be  equal  to  or  greater  than  the  torque  provided  by  the  bending 

moment,  i.e.f  Tb  =  (£'/b)Mb  or  Tbe/Tb  ^  1,  and  from  (C-2)  and  (C-12): 

Tbe/Tb  =  [3kiE/ob(l+v)(b/Jl')]bd2[(d/LT)4>s  +  (d/*’)$p]  l>  1.0,  or 

n  =  (bA')[3k1(E/ob)/(ltv)]((d/LT)4s  +  Cd/fi.«)4>F]  i  1-0  (C-13) 

where  ab  is  the  bending  stress  at  failure  according  to  (C-l). 

Note  that  Tb/Mj,  =  b/s,'  and  for  four-point  loading  £  *  =  a;  for  three-point 
loading  £ '  =  L/2,  thus  (C-10)  becomes: 

Xax  =  (°b/2)  (d+v)  +  (l/3k2)[(l/(«,'/b))2  +  9(l-v)2k22]^|  (C-14) 

with  n  =  1.0. 

Case  II 

If,  however,  n  <  1.0  then  bottoming  occurs  prior  to  or  at  failure  and  the 
following  analysis  is  applicable: 

Equations  (C-4)  and  (C-5)  are  still  appropriate  but  the  shear  stress  is 
t  =  Tbe/k2bd2.  (C-15) 

Substitution  of  x  from  the  above  and  ab  from  (C-l)  into  (C-5)  gives 
tan28*  =  (Tbe/Mb)/3(l-v)k2, 
but  from  (C-13)  Tbg  =  nTb  and  thus: 

tan26*  =  (nTbe/Mb)/3(l-v)k2.  (C-16) 

As  in  Case  I,  using  a  like  procedure  we  determine  anmax  to  ^e: 

°nmax  =  (a*/2)  +  (l/3k2)  [n2(Tb/Mb) 2  +  9(l-v)2k22]^  (C-17) 

with 

n  =  [3k1(E/crb)/(l+v)][(d/Lr)4.s  +  (d/t')*F]  <  1.0-  (C-18) 

Equation  (C-17)  is  applicable  to  both  systems  since  (Tb/Mb)  =  b/ £ 1  =  b/a  for  the 
three-point  beam  bending  system,  and  (Tb/Mb)  =  £'/b  =  2b/L  for  the  four-point 
system.  Thus  (C-17)  becomes 

°nmax  =  (0b/2)  {^1+v)  +  f1/3^)  f  (nb/^)2  ♦  9(l-v)2k22]54  (C-19) 


jfc.  amitafettB  2m 
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Notice  when  n  =*  1,  (C-19)  reduces  to  (C-14) . 

Finally,  the  percent  error  is  defined  as: 

F  =  t^b^nmax)/^^^100*  CC‘20) 

Errors  were  calculated  in  accordance  with  (C-20) ,  with  v  =  0.25  for  Case  I  (n  = 
0.20,  0.40,  0.60,  and  0.80)  and  Case  II  (n  =  1.0).  These  results  are  given  in 
tabular  form  as  a  function  of  b/d  and  l'/b  in  Tables  10a  to  lOe  in  the  text. 


APPENDIX  D._  WEDGING  STRESSES 

We  allow  the  stress  in  the  x  direction  in  Figures  1  and  2  to  be 

ax  =  ob  +  (2P/bd)&r,  (D-l) 

where  %  is  the  bending  stress,  i.e.,  crb  =  (6Mx/bd2) ,  and  2P/bd  is  the  local 
stress,  i.e.,  the  so-called  wedging  stress,  and  Bt  is  a  numerical  factor  dependent 
on  the  normalized  distance  x'/(d/2)  on  either  side  of  the  applied  load  point.2 

For  convenience  the  value  of  By  at  the  tensile  side  of  the  beam  as  a  function 
of  x'/(d/2)  is  given  in  Table  D-l. 

The  percent  error  is  defined  as: 

e-  [(ab-ax)/ox]100,  CD-2) 

and  substitution  of  (D-l)  into  the  above  equation  gives 

e=  ^-BT/[(ab/2P/bd)+BT]  100.  (D-3) 

Since  ab  =  6Mx/bd2,  then  (D-3)  becomes: 

7  =  ^-eT/(3Mx/Pd+Bi)}  100.  (D-4) 

For  a  four-point  loaded  beam  the  bending  moment  is  constant,  i.e.,  Mx  =  Pa, 
and  thus  equation  (D-4)  becomes: 

e  =  {-8T/[(3a/d)+BT]}l00.  (D-5) 

From  Table  D-l  and  Equation  D-5  above,  it  is  seen  that  the  error  is  dependent 
on  or  the  fracture  location,  which  is  the  normalized  distance  x'/(d/2).  These 
errors  have  been  computed  for  the  four-point  loaded  beam  and  presented  in  Table  11a. 

For  the  three-point  loaded  beam,  (D-4)  is  still  applicable,  but  recalling  that 
Mx  =  P/2[ (L/2) -x ']  and  substituting  Mx  into  (D-4)  gives: 

e  ={-gT/(i(L/d)-2(x7d)*BT)}l00.  (D-6) 
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Again,  as  can  be  seen  by  (D-6) ,  the  percent  error  is  dependent  upon  and 
the  normalized  fracture  location.  These  percent  errors  are  given  in  Table  lib. 


Table  D-l.  WEDGING  PARAMETER  eT 


±x  '/d 

S*T 

0 

-0.1332 

0.125 

+0.0137 

0.250 

+0.0868 

0.375 

+0.0640 

0.500 

+0.0421 

0.750 

+0.0220 

1.000 

+0.0095 

1.500 

+0.00075 

APPENDIX  E.  CONTACT  POINT  TANGENCY  SHIFT 

Consider  the  four-point  loaded  beam  shown  in  Figure  4.  The  original  span  "a 
is  seen  to  decrease  by  the  amount  (hi +h2  )  due  to  rolling  or  slipping  of  the  beam 
on  its  support  and  load  points.  The  beam  fulfills  the  condition:* 

dy2/dx2  =  MX/EI.  (E-l) 

The  moments  are  defined  as: 

Mx  =  P(x-h|),  0  <  x  <  (a-h2) 

Mx  =  P[a-(h1+h2) ] ,  (a-h2)  <  x  <  L-(a-h2). 

0  <  x  £  (a-hj) ,  and  (E-2) 

(a-h2)  <  x  <  L-(a-h2)  .  (E-3) 

=  P[a-(hj+h2) ] (a-h2)  +  C2,  or 

C2-Cx  =  -(P/2) (a-h2)2.  (E -4 ) 

Note  also  that  when  x  =  L/2  and  dy/dx  =  0  in  (E-3), 

C2  =  -P[a-(h1+h2)]L/2.  (E5) 

After  substitution  of  C2  into  (E-4)  we  obtain: 

Cx  =  (P/2)  ^(a-h2)2  -  L[a-(hx+h2)]j.  .  (E-6) 

•It  is  realized  that  a  change  in  span  dimensions  must  give  rise  to  horizontal  reactions;  but  it  is  assumed  that  this  effect  is  small, 
and  is  thus  ignored. 


The  slope  equations  are 

El (dy/dx)  =  P[(x2/2)-hjx]  +  Cj, 

EI(dy/dx)  =  P[a- (hi+h2) ]x  +  C2, 

Now  when  x  =  a-h2, 

Cx  +  P  {[(a-h2)2/2]  -  (a-h2)hx| 
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Substitution  of  these  constants  into  the  appropriate  slope  equations  gives: 

El (dy/dx)  =  P[(x2/2)-hlX]  +  (P/2)  {(a-h2)2  -  L[a-(hi+h2)] J  (E-7) 

with  0  <  x  <  (a-h2),  and 

El (dy/dx)  =  P[a-(h1+h2)](x-(L/2)],  (E-8) 

with  (a-h2)  <  x  <  L-(a-h2). 

Now  when  x  =  hi  from  the  geometry  dy/dx  =  -hx/Af-hf  »  -hi/Ri  and  when 
x  =  a-h2,  dy/dx  =  -ha/i^F-hF  *  -h2/R2*  The  above  relationships  are  used  with  (E-7) 
and  (E-8)  and  we  obtain: 


(Rl/d)  = 


(h^a)  (E/ob) 


(hi/a)2  -  (l-h2/a) 2  +  (L/a)[l-(hi/a  +  h2/a)] 


(E — 9) 


and 


where 


1-  (h2/a)  =  (1/2)  {[(L/2a)  +  (hi/a)  +  A2] 

-  «'[(L/2a)  +  (hi/a)  +  A2]2  -  4 [(hj/a)  (L/2a)^  A2]| 
A2  =  (E/2ob)/(R2/d). 


Note  that  for  the  three-point  loading  case  h2  =  0,  L/a  *  2,  and  P 
(E-9)  reduces  to: 

Ri/d  =  [2(hi/L)(E/ab)]/[(2hi/L)-l]2, 
and  the  region  of  validity  of  (E-10)  vanishes. 

The  percent  error  is  defined  as: 

e  -  [(%-ax)/ax]100  =  [(Mb-Mx)/Mx]100,  or 

i  =  |[(hi/a)  +  (h2/a)]/[l  -  (hj/a)  -  (h2/a)]} 

for  four-point  loading  and 

i  =-j(2h1/L)/[l  -  (2hi/L)]jl00 

for  three -point  loading. 


(E-10) 

P/2  and 
(E —11) 


100 


(E-12) 


(E- 13) 


Calculations  were  performed  by  the  following  procedure.  It  was  assumed  that 
E/ob  =  1  x  103,  and  thus  A2  =  (1  x  103)/(p2/d);  then  for  the  1/3-  and  1/4-four- 
point  loading  case  a/L  =  1/3  and  1/4,  numerical  values  were  assigned  to  p2/d  and 
hi/a  and  then  h2/a  was  determined  from  (E-10).  Those  numerical  values  of  hi/a 
and  corresponding  h2/a  were  substituted  into  (E-9)  to  determine  pi/d.  This  same 
procedure  was  used  for  the  three-point  loading  case  with  h2  =  0  and  L/a  =  2.  Once 
the  parameters  hi/a,  h2/a,  p2/d  and  pi/d  are  known,  percent  errors  according  to  (E-12) 


for  the  four-point  loading  case  and  (E-13)  for  the  three-point  loading  case  can  be 
determined.  Such  errors  are  given  in  Tables  12a-c. 


APPENDIX  F.  ERROR  DUE  TO  NEGLECTING  CHANGE  IN  MOMENT  OF  INERTIA 
CAUSED  BY  CORNER  RADII  OR  CHAMFERS 

Corner  Radius 

Consider  Figure  Sa,  which  shows  the  cross  section  of  a  rectangular  beam  with 
corner  radii  r.  The  true  moment  of  inertia  (Ixx)r  about  the  centroidal  or  neutral 
axis  x-x  is: 

(Xxx)r  =  b(d-2r)3/12  +  (b-2r)r3/6  +  (1/2) (b-2r) (d-r) 2r 

+  2r‘t(ir/9  -  8/9tt)  +  irr2[d/2  -  r(l  -  4/3ir)]2.  (F-l) 

Most  investigators,  however,  will  neglect  the  loss  of  inertia  when  calculating  the 
bending  stress  due  to  corner  radii  and  assume  that  1^  =  bd3/12,  and  if  we  define 
the  percent  error  as  { [I^-(Ixx)r]/(Ixx)r}100,  we  then  obtain: 

1  -  {(l-2r/d)3  +  2[l-2(r/d) (d/b) ] (r/d) 3 
+  6[l-2(r/d)  (d/b)  ]  (1 -r/d) 2 (r/d)  +  24 (d/b)  (r/d)  4(*/9-8/9tt) 

+  12rr(d/b)  (r/d)2[l/2-(r/d)  (l-4/3n)]2| 

(l-2r/d) 3  +  2 [1-2 (r/d) (d/b) ] (r/d) 3  ♦  6 [1 -2 (r/d) (d/b) ] (1-r/d) 2 (r/d) 

+  24  (d/b)  (r/ d) 14  (ir/9-8/9ir)  +  12*  (d/b)  (r/d)  2  [1/2-r/d  (l-4/3ir)  ] 2  . 

(F-2) 

45'  Chamfer 

Now  consider  Figure  5b,  which  shows  a  rectangular  beam  with  45°  corner  cham¬ 
fers  c.  The  true  moment  of  Inertia  (Ixx)c  about  the  x-x  axis  is: 

(Ixx)c  =  (bd3/12)  -  (c2/9)  [c2  ♦  (1/2) (3d-2c)2  j  .  (F-3) 

Defining  the  error  in  the  same  manner  as  above  results  in: 
r (4/3) (d/b) (c/d) 2 [(c/d) 9  +  (1/2) (3-2c/ d) 2 ] 

'  ~  -  (4/3)  (d/b)  (c/d) 2  [(c/d) 2  +  (1/2)  (3- 2c /d)2] 

The  errors  were  calculated  for  various  values  of  d/b  as  a  function  of  r/d 
from  (F-2)  and  c/d  from  (F-4) .  These  results  are  shown  in  Tables  13  and  14. 


(F-4) 
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